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ABSTRACT 

We present modeling and interpretation of the continuum and emission lines for a 
sample of 51 unobscured Type 1 active galactic nuclei (AGN). All of these AGNs have 
high quality spectra from both XMM- Newton and Sloan Digital Sky Survey (SDSS). 
We extend the wavelength coverage where possible by adding simultaneous UV data 
from the OM onboard XMM-Newton. Our sample is selected based on low reddening 
in the optical and low gas columns implied by their X-ray spectra, except for one case, 
the BAL-quasar PG 1004+130. They also lack clear signatures for the presence of a 
warm absorber. Therefore the observed characteristics of this sample are likely to be 
directly related to the intrinsic properties of the central engine. 

To determine the intrinsic optical continuum we subtract the Balmer continuum 
and all major emission lines (including Fell). We also consider possible effects of 
contamination from the host galaxy. The resulting continuum is then used to derive 
the properties of the underlying accretion disc. We constrain the black hole masses 
from spectral fits of the Balmer emission lines and determine the best fit value from 
the modeling of broadband spectral energy distributions (SED). In addition to the 
disc component, many of these SEDs also exhibit a strong soft X-ray excess, plus a 
power law extending to higher X-ray energies. We fit these SEDs by applying a new 
broadband SED model which comprises the accretion disc emission, low temperature 
optically thick C omptonisation an d a hard X-ray tail by introducing the concept of 
a corona radius ()Done et al.ll201lT) . We find that in order to fit the data, the model 
often requires an additional long wavelength optical continuum component, whose 
origin is discussed in this paper. We also find that the Photo-recombination edge of 
Balmer continuum shifts and broadens beyond the standard limit of 3646A, implying 
an electron number density which is far higher than that in the broad line region 
clouds. 

Our results indicate that the Narrow Line Seyfert Is in this sample tend to have 
lower black hole masses, higher Eddington ratios, softer 2-10 keV band spectra, lower 
2-10 keV luminosities and higher a OXl compared with typical broad line Seyfert Is 
(BLS1), although their bolometric luminosities are similar. We illustrate these differ- 
ences in properties by forming an average SED for three subsamples, based on the 
FWHM velocity width of the H/3 emission line. 

Key words: accretion, broadband SED modeling, active-galaxies: nuclei 



1 INTRODUCTION 

The spectral energy distribution (SED) of AGN has been 
modeled for several decades. Initial studies foc used on the 
infra r ed, optical and ultraviolet continuum (e. g. IWills et al.l 



1 19851 ; ICanalizo fc Stockton! l200ll ; lLacv et all [20071 ). With 



E-mail: chichuan.jin@durham.ac.uk 



the inclusion of X-ray data, it was possible to define the 
continuum on both sides of the ultraviolet/X-ray gap (im- 
posed by galactic photoelectric absorption ), and so constrain 
the properties o f the accretion disc (e.g. IWard et al.l 1 19871 : 
lElvis et al.lll994l y Refinements to modeling the optical/UV 
continuum include subtraction of the complex blended fea- 
tures arising from permitted iron emission, the so-called 
small blue-bump from the Balmer continuum, and contami- 
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nation across the entire spectrum from a stellar component 
l|Maoz et al.lll993l ; iBoisson et aI]|200Ch 

The observed spectral differences between various types 
of AGN are not only due to selective absorption and 
orientation effects, as i mplied by the sim plest version of 
AGN unification model i|Antonuccilll993l ), but also result 
from a wide range in basic physical pa rameters, such as 
black hole mass and accretion rate (e.g. Boroson fc Greer] 
19921; IBoller, Brandt fc Finklll996l ; iDone fc Gierliriskill2005l ; 
Zhou et alj 20061 ). To better understand the accretion pro- 
cesses occurring close to the super massive black hole 
(SMBH), we construct broadband SEDs. Galactic dust red- 
dening, together with the intrinsic reddening of the AGN it- 
self, attenuates the optical/UV band emission. Furthermore, 
Photoelectric absorption from gas modifies the lower energy 
X-ray continuum. But these factors can be quantified and 
corrected. Thereby we can recover the intrinsic SED, except 
for the unobservable far-UV region. If we have reliable data 
on both sides of the energy gap between the UV and soft 
X-ray, we can apply a multi-component model which spans 
across it. 



1.1 Previous Work 

Many multi-wavelength studies h ave been carried out previ- 
ously. IPuchnarewicz et al. 1 19921) s t udied the optical prop- 
erties of 53 AGNs in ICordova et ail (l992)'s sample with ul- 
tra soft X-ray excesses, and found that they tend to have 
narrower permitted lin es than optically selected sa mples. 
Supporting this finding. IBoller. Brandt fc Finkl l|l996l ) stud- 
ied ROSAT selected AGN with extremely soft X-ray spec- 
tra, and found that they tend to be Narrow-Line Seyfert Is 
(NLSls). Correspondingly they found that o ptically selected 
NLSl s often have large soft X-ray excesses. IWalter fc Finkl 
|l993t ) combined soft X-ray and optical data for 58 Seyfert 
Is, and showed that their broadband SED have a bump 
from UV to soft X-r ays, which is now refer ed to as the 
big bl ue bump ('BBB). lGrupe et al.l (|l998h and lGrupe et all 
ljl999h used a sample of 76 bright soft X-ray selected Seyferts 
with infrared data, optical spectra and soft X-ray spectra. 
Their results reinforced the connection between the opti- 
cal and soft X-ray spectra, and confirm ed the exist e nce o f 
strong BBB emission in these objects. lElvis et al l (|l994l ) 
studied 47 quasars in a UV-soft X-ray sample, and derived 
the mean SEDs for radio-loud and radio-quiet sources. Re- 
cently, more detailed spectral models have been applied to 
broadband SEDs including simultaneous optical/UV and 
X-ray observ ations which avoid po t entia l problems caused 
by variability. I Vasudevan fc Fabianl ()2007i ) (hereafter VF07) 
combined a disc and broken powerlaw model to fit opti- 
cal, far UV and X-ray data for 54 AGN. They found a 
well-defined relationship between the har d X-ray bolomet- 
ric co rrection and the Eddington ratio. iBrocksopp et al] 
(2006) analysed the data from XMM- Newton's simultane- 
ous EPIC (X-ray) and OM (optical/UV) observations for 
22 Palomar Green (PG) quasars. Another sample consist- 
ing of 21 NLSls and 13 broad line AGNs was also defined 
using simultan eous data from XM M-Newton's EPIC and 
OM monitor jCrummv et al.ll2006l ). The SEDs of this sam- 
ple were then fitted using various broadband SED models 
such as disc plus powerlaw m odel, disc reflection model and 
disc wind absorption model (Mi ddleton. Done fc Gierlinskil 



120071 ). IVasudevan fc Fabianl (|2009h derived SEDs using 
XMM-Newton's simultaneous X-ra y and optical/UV obse r- 
vations for 29 AGNs selected from iPeterson et all (|2004l )'s 
reverberation mapped sample. The well constrained black 
hole masses available for this sample enabled them to fit 
a better constrained accretion disc model, combined with 
a powerlaw, to the source's broadband SEDs. Hence they 
derived more reliable Eddington ratios. 



1.2 Our AGN Sample 

In this paper we define an X-ray/optically selected sample 
of 51 AGN, all of which have low reddening (so excluding 
Seyfert 2s and 1.9/1. 8s), to construct SEDs ranging from 
about 0.9 microns to 10 keV. We also apply corrections for 
the permitted iron features, the Balmer continuum and stel- 
lar contribution, in order to model the non-stellar continuum 
free from emission line effects. Included in this sample are a 
number of NLSls, a subclass of AGN whose permitted line 
widths are comparable to those of forbidden lines. Their 
[OIII]A5007/H/3 ratio is also lo wer than the typical value of 
broad line Seyfert Is (BL Sls) (|Shuder fc Osterbrockl Il98ll ; 
lOsterbrock fc Poggel Il985l ) . For consistency with previous 
work, we classify AGNs in our sample as NLSls if they have 
ratios of [OIII1A 5007/H/3 < 3 and FWHM H/3 < 2000 km/s 
l|Goodrichll 19891 ). We identify 10~12 NLSls in our samplqj. 

All objects in our sample have high quality optical spec- 
tra taken from the Sloan Digital Sky Survey (SDSS) DR7, 
X-ray spectra from the XMM-Newton EPIC cameras, and 
in some cases simultaneous optical/UV photometric data 
points from the XMM-Newton OM monitor. Combining 
these data reduces the impact of intrinsic variability and 
provides a good estimate of the spectral shape in the opti- 
cal, near UV and X-ray regions. In addition, by analyzing 
the SDSS spectra, we can derive the parameters of the prin- 
cipal optical emission lines and underlying continuum. An 
important result from reverberation mapping study is the 
correlation between black hole mass, mo nochromatic lumi- 
nosity at 5100 A and R/3 FWHM (e.g . iKaspi et all 120001 ; 
IWoo fc Urrvll2002l ; IPeterson et al.ll2004r i. We measure these 
quantities from the SDSS spectra, and then estimate black 
hole masses using this correlation. 

Compared with previous work, a significant improve- 
ment of our study is that we employ a new broadband SED 
model which combines disc emission, Comptonisation and a 
high energy powerlaw component in the context of an en- 
erget ically self-consiste nt model for the accretion disc emis- 
sion (|Done et alj|201ll . also see Section f5 . 2 f) . By fitting this 
model to our data, we can reproduce the whole broadband 
SED from the optical to X-ray. From this detailed SED fit- 
ting, we derive a number of interesting AGN properties such 
as: the bolometric luminosity, Eddington ratio, hard X-ray 
slope, and the hard X-ray bolometric correction. Combin- 
ing all the broadband SED parameters with the optical pa- 
rameters, we can provide further evidence for many previ- 



1 Although 2XMM J112328.0+052823 and IE 1346+26.7 have 
H/3 FWHMs of 2000 km s~ 1 , 2050 km s~ 1 respectively, they both 
have Ha FWHM of 1700 km s" 1 , and also share other NLSl's 
spectral characteristics. Thus they could both potentially be clas- 
sified as NLSls, making a total of 12. 
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Table 1. The Scyfcrt 1 Galaxy Sample Set 
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ID 
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Redshift 
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Obs Date 
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Obs Date 


Counts c 
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2000-09-25 


2056 


3 
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2001-10-23 
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2XMMi/DR7 
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2001-04-26 
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2002-12-10 
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5 


2XMMi/DR7 


0.358 


J080608. 0+244421 


2001-10-26 


52705-1265-410 


2003-03-07 


2912 


6 
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J081422. 1+514839 


2003-04-27 


53297-1781-220 


2004-10-19 
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7 
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J092246. 9+512037 
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9 
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2004-11-15 
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2001-12-09 
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10 


PG 0947+396 


0.206 


J095048. 3+392650 


2001-11-03 


52765-1277-332 


2003-05-06 


58555 


11 
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J100025. 2+015852 


2003-12-10 


52235-0501-277 


2001-11-22 
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12 
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2004-04-20 
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2003-02-02 
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13 


PG 1004+130 


0.241 


J100726. 0+124856 
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53055-1744-630 


2004-02-20 
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14 


RBS 0875 


0.178 
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15 
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16 
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0.160 
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J 105143. 8+335927 
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17 


1RXS J111007 
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Jll 1006.8+612522* 


2006-11-25 


52286-0774-600 


2002-01-12 


6147 


18 


PG 1115+407 


0.155 


Jl 11830. 2+402554 


2002-05-17 


53084-1440-204 


2004-03-20 


64601 


19 


2XMM1/DR7 


0.101 


J 112328.0+052823 


2001-12-15 


52376-0836-453 


2002-04-12 


10098 


20 


RX J1140. 1+0307 


0.081 


J114008. 7+030710 


2005-12-03 


51994-0514-331 


2001-03-26 


35616 


21 
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J 120442. 1+275412 


2002-05-30 
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2006-03-25 
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22 
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2001-06-15 
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2005-04-06 
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23 
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J121930. 9+064334 
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2004-04-26 
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24 
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J122018.4+064120 


2002-07-05 
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25 
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J 123054. 1+1 10011 


2005-12-17 


52731-1232-417 


2003-04-02 


1 c f o n o 

165823 


26 


2XMMi/DR7 


0.304 


J123126. 4+1051 11 


2005-12-17 


52731-1232-452 


2003-04-02 


8816 


27 


MRK 0771 


0.064 


J123203. 6+200929 


2005-07-09 


54481-2613-342 


2008-01-15 


40705 


28 


RX J1233. 9+0747 


0.371 


J123356. 1+074755 


2004-06-05 


53474-1628-394 


2005-04-14 


6041 


29 


RX J1236. 0+2641 


0.209 


J123604. 0+264135* 


2006-06-24 


53729-2236-255 


2005-12-25 


17744 


30 


PG 1244+026 


0.048 


J124635. 3+022209 


2001-06-17 


52024-0522-173 


2001-04-25 


8509 


31 


2XMMi/DR7 
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2000-06-21 
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2006-03-26 
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32 
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2004-06-06 
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2005-05-09 
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33 
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J132101. 4+340658 


2001-01-09 


53851-2023-044 


2006-04-26 


4425 


34 


1RXS J 132447 
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J132447.6+032431 


2004-01-25 


52342-0527-329 


2002-03-09 


6305 


35 
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J134113. 9-005314 


2005-06-28 


51671-0299-133 


2000-05-07 


18007 


36 
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0.059 


J134834. 9+263109 


2000-06-26 


53848-2114-247 


2006-04-23 


71985 


37 


PG 1352+183 


0.151 


J135435. 6+180518 


2002-07-20 


54508-2756-228 


2008-02-12 


36171 


38 


MRK 0464 


0.050 


J135553. 4+383428 


2002-12-10 


53460-2014-616 


2005-03-31 


13974 


39 


1RXS J135724 


0.106 


J135724.5+652506 


2005-04-04 


51989-0497-014 


2001-03-21 


12081 


40 


PG 1415+451 


0.114 


J141700.7+445606 


2002-12-08 


52728-1287-296 


2003-03-30 
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41 


PG 1427+480 


0.221 


J142943.0+474726 


2002-05-31 
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2005-04-01 


70995 


42 


NGC 5683 


0.037 


J143452.4+483943 


2002-12-09 


52733-1047-300 


2003-04-04 
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43 
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0.208 
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2005-02-11 
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2005-05-04 
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41 
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0.065 
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2003-02-08 


54208-2142-637 


2007-04-18 


134532 


45 
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2002-08-25 


53083-1353-580 


2004-03-14 
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46 
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2001-08-21 


54563-1835-054 


2008-04-07 


10952 


47 
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0.090 


J155829.4+271715 


2002-09-10 


52817-1391-093 


2003-06-27 


6995 


48 


MRK 0493 


0.031 


J155909.6+350147 


2003-01-16 


53141-1417-078 


2004-05-14 


124115 


49 


II Zw 177 


0.081 


J221918.5+120753 


2001-06-07 


52221-0736-049 


2001-11-08 


36056 


50 


PG 2233+134 


0.326 


J223607.6+134355 


2003-05-28 


52520-0739-388 


2002-09-03 


7853 


51 


MRK 0926 


0.047 


J230443.3-084111 


2000-12-01 


52258-0725-510 


2001-12-15 


59513 



a for some targets without well-known names, we simply use '2XMMi/DR7'; 

b the full name should be '2XMM J...', but for those targets with * symbol, their full names should be '2XMMi J...'; 

c the total counts in all three EPIC monitors, namely pn, MOS1 and MOS2, and there are at least 2000 counts in at least one 

of these three monitors; 
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ously suggested correlations, including all the correlations 
between optical and X-ray claimed in previous work, plus 
many others such as the H/3 FWHM versus X-ray slope, 
black hole mass versus Eddington ratio, Fell luminosity ver- 
sus [OIIIJA5007 emission line luminosity and the high excita- 
tion lines (e. g. [FeVII1A6087, [FeX]A6374) versus their ioniz- 
ing flux (e.g. Boroson fe Green 19921; Boiler. Brandt fc Fink 



19961 ; IGrupe et all 119981; IGrupe et al.1 Il999l ; ISulentic et al 
200ol ; iMullanev et al.ll2009h . 

This paper is organized as follows. Section 2 describes 
the sample selection and data analysis procedures. The de- 
tailed spectral fitting methods and results including Balmer 
line fitting, optical spectral fitting and broadband SED fit- 
ting are each discussed in sections 3, 4 and 5, separately. We 
present the statistical properties of our sample in section 6. 
The summary and conclusions are given in section 7. A flat 
universe model with Hubble constant of Ho = 72 km s _1 
Mpc -1 , Qm = 0.27 and J7a = 0.73 is adopted. In another 
paper, we will present our analysis of correlations between 
selected optical/UV emission features and the SED compo- 
nents, and discuss their physical implications (Jin et al. in 
prep., hereafter Paper II). 



2 SAMPLE SELECTION AND DATA 
ASSEMBLY 

To identify a sample of Type 1 AGNs having both high 
quality X-ray and optical spectra, we performed a cross- 
correlation between 2XMMi catalog and SDSS DR7 catalog. 
We filtered the resulting large sample as described below. 
Our final sample consists of 51 Type 1 AGNs including 12 
NLSls, all with high quality optical and X-ray spectra and 
low reddening/ absorption, and with H/3 line widths ranging 
from 600 fcms" 1 up to 13000 ferns -1 . All the sources are 
listed in Table [t] 



2.1 The Cross-correlation of 2XMMi & SDSS 
DR7 

The first step was to cross-correlate between 2XMMi and 
SDSS DR7 catalogs. The 2XMMi catalog contains 4117 XMM- 
Newton EPIC camera observations obtained between 03-02- 
2000 and 28-03-2008, and covering a sky area of ~ 420 deg 2 . 
The SDSS DR7 is the seventh data release of the Sloan sky 
survey. The SDSS spectroscopic data has sky coverage of ~ 
8200 deg 2 , with spectra from 3800 A to 9200 A, and spectral 
resolution between 1800 and 2200. 
Our cross-correlation consisted of three steps: 

1. We first searched for all XMM/SDSS position pairs that 
lay within 20" of each other, resulting in 5341 such cases. 

2. For these 5341 unique X-ray sources, we imposed two fur- 
ther selection criteria: that source positions be separated by 
less than 3", or that sources be separated by no more than 
3 x the XMM-Newton position uncertainty and no more 
than 7". This filtering resulted in 3491 unique X-ray sources. 
The 3" separation is chosen because we want to include all 
possible XMM/SDSS pairs during these early filtering steps. 
From the 2XMMi and SDSS DR7 cross-correlation, there are 
114 XMM/SDSS pairs whose separations are less than 3", 
but are still nevertheless greater than 3 x the XMM position 



uncertainty. We included all of these pairs. The 7" separa- 
tion upper limit mitigates spurious matches, especially for 
fainter objects and/or those located far off-axis. 
3. We selected only objects classified as extragalactic, giving 
a total of 3342 for further analysis. 



2.2 Selection of Seyfert 1 with High Quality 
Spectra 

Within these 3342 unique X-ray sources which satisfied all 
the above criteria, we applied further filtering to select only 
Type 1 AGNs having both high quality optical and X-ray 
spectra. The five steps in the filtering were as follows: 

1. In order to obtain black hole mass estimates and also as 
reddening indicators, we require H/3 and Ha emission lines 
to be measurable. So we only selected sources with H/3 in 
emission (as indicated by the SDSS H/3 line models with at 
least 3cr significance and EW > 0) and redshift z < 0.4. 
This selection resulted in 802 unique X-ray sources, and 888 
XMM/SDSS pairs (since some X-ray objects were matched 
with more than one SDSS spectrum). 

2. Then we searched for the Type 1 AGNs (including sub- 
types 1.0, 1.5, 1.8 and 1.9) which have a minimum of 2000 
counts in at least one of the three EPIC cameras. Our search 
retrieved 96 such broad line AGNs. We then inspected each 
of these XMM/SDSS pairs, to confirm that all the matches 
were indeed genuine. 

3. From inspection of the SDSS spectra, we excluded 22 
sources whose blueward part of the H/3 line showed strong 
reddening or low S/N, which would distort the H/3 line pro- 
file. We also excluded one object, RBS 0992, because its 
SDSS spectrum did not show an H/3 line, due to a bad data 
gap. We ensured that the remaining 73 objects all had good 
H/3 line profiles. 

4. As a simple method to assess the spectral quality of the X- 
ray data, we used wabs*powerlaw model in xspecll.3.2 to 
fit the rest-frame 2-10 keV X-ray spectra of all 73 objects. 
The error command was used to estimate the 90% confi- 
dence region for the photon index parameter. Based on the 
results, 16 objects with photon index uncertainties greater 
than 0.5 were thereby excluded, leaving 57 Type 1 AGNs 
with relatively well constrained 2-10 keV spectra. 

5. By examining the 0.2-10 keV X-ray spectra, we ex- 
cluded another 6 objects (i.e. IRAS F09159+2129, IRAS 
F12397+3333, PG 1114+445, PG 1307+085, PG 1309+355 
and PG 1425+267) whose spectral shapes all showed 
clear evidence of an absorption edge at ~0.7 keV (pos- 
sibly originating fro m combined Fe I L-Shell and O V II 
K-Shell absorptions (|Lee et al.1 lioOll ; iTurner et al.ll2004h ). 
This is a typical spectral signature of a warm absorber 
(e.g.)Nandra fc Pounds! 19941 ; [Crenshaw. Kraemer fc George! 
120031 ). By removing such objects with complex X-ray spec- 
tra, our broadband SED fitting is simplified. Our final sam- 
ple contains 51 Type 1 AGNs. 



2.3 Characteristics of the Sample 

The sample selection procedure described above ensures that 
every source in our AGN sample has both high quality op- 
tical and X-ray spectra. In addition, a large fraction of the 
sample have simultaneous optical/UV photometric points 
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Figure 1. The aperture effect correction results for 17 extended sources in the sample. The point like source RBS 0769 (the last figure 
marked by **) is also shown for comparison. We over-plot OM data points on to the SDSS spectrum. Red OM points are data obtained 
directly from the OM PPS files. Blue OM points are the corresponding data after applying a smaller 6" aperture to all OM filters, and 
applying appropriate OM corrections to the flux eg. deadtime correction, coincidence loss correction and OM time sensitivity degradation 
correction. 
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from the OM monitor. Such high quality data enables ac- 
curate spectral fitting. In the optical band our sample is 
selected to have low reddening, since if present this would 
significantly modify the intrinsic continuum as well as the 
optical emission lines. This requirement reduces the com- 
plexity and uncertainty in our modeling of the intrinsic con- 
tinuum, and also increases the overall quality of H/3 and 
Ha line profiles useful for estimating the black hole masses. 
Furthermore, low reddening is essential in the UV band. The 
inclusion of OM-UV photometric data observed simultane- 
ously with the X-ray spectra provides a reliable link between 
these bands. This helps to reduce fitting uncertainty of the 
SED resulting from optical and X-ray variability. Besides, 
all sources are well constrained in the 2-10 keV band, which 
is directly associated with the compact emitting region of 
the AGN. Our exclusion of objects with evidence of a warm 
absorber means that the 2-10 keV spectral index is likely to 
be intrinsic rather than hardened by absorption in the soft 
X-ray region. 

In summary, compared with previous AGN samples 
used for broadband SED modelling, the spectrally 'cleaner' 
nature of our sample should make the reconstructed broad- 
band SEDs more reliable. Consequently, the parameters de- 
rived from the broadband spectral fitting should be more 
accurate. This may reveal new and potentially important 
broadband correlations, which we will discuss in detail in 
paper II. 

2.4 Additional Data 

The 51 Type 1 AGNs all have SDSS survey- quality spectra 
(flagged as "sciencePrimary" in SDSS catalog), including 3 
objects that have multiple SDSS spectra (i.e. NVSS J030639, 
1RXS J111007 and Mrkl018). In such cases we adopt the 
SDSS spectrum which connects most smoothly with the OM 
data. 

For each object, we used all available EPIC X-ray spec- 
tra (i.e. pn, MOS1 and MOS2) for the broadband SED mod- 
eling, unless the spectrum had few counts and low S/N. We 
also searched through the XMM-OM SUSS catalog for all data 
in the OM bands (i.e. V, B, U, UVW2, UVM2 and UVW1), 
which are observed simultaneously with the corresponding 
EPIC spectrum. Of our 51 sources, we have 14 sources with 
SDSS optical spectra and XMM EPIC X-ray spectra, and 
37 sources which in addition to this also have XMM-OM 
photometry. 



2.5 OM Data Corrections and Aperture Effects 

In the procedure of combining the SDSS spectra and OM 
data points, we identified that in some objects there is a 
clear discrepancy between these two data sets. The OM 
points often appear higher on the spectral plots (brigher) 
than is consistent from a smooth extrapolation of the SDSS 
spectral shape. In fewer cases this discrepancy appears in 
the opposite sense, with the OM points apparently too low 
(fainter), see Figure [T] for some examples). This discrepancy 
may arise for several reasons, including a simple aperture 
effect. Compared to 3" diameter for the SDSS spectroscopy 
fibres, the OM monitor has a much larger aperture, i.e. 12" 
and 35" diameter for the OM optical and OM UV filters 



respectively l|Antonio Talavera.OMCal Team! 120091 ). If the 
host galaxy is sufficiently extended, e.g. in the case of RE 
J 1034+396, the larger aperture of the OM would include 
more host galaxy emission than that in the SDSS spectrum 
(see also section [5.3.11 for other possible reasons to account 
for this discrepancy). To investigate the aperture issue in 
more detail, we performed the following tests: 

(1) We examined the combined SDSS and OM data plots, 
searching for those objects with excess OM flux compared 
with that expected from the extrapolated SDSS spectrum. 
We identified 27 such cases out of the 51 sources; 

(2) Within this sample of 27 sources, we checked the cata- 
log flag for an extended source in each OM filter. We noted 
those flagged as an extended source in at least one OM fil- 
ter. This yielded 13 sources out of the 27. 

(3) We also extracted the SDSS CCD images for all 51 ob- 
jects and visually checked whether they appeared extended. 
As a result, we included another 4 objects for which their 
SDSS CCD images show that their host galaxy is extended 
beyond the 3" diameter of the SDSS aperture. Either they 
were not flagged as extended sources in any OM filter, or 
they did not have any OM optical data. For these 17 objects, 
an aperture effect could at least be partially responsible for 
an excess flux in the OM data. 

(4) For these 17 objects we downloaded all available OM 
image files. In each OM image, we applied a 6" diameter 
aperture from which to extract the flux. We used the same 
sized aperture placed on a blank region of sky close to the ob- 
ject, to estimate the background. The quoted PSF FWHM 
of the OM for the different filters are: V(1.35"), B(1.39") , 
U(1.55"), UVW1(2.0"), UVM2(1.8"), UVW2(1.98"). Thus 
in all cases 6" is at least 3xPSF FWHM. So this aperture 
includes effectively all optical flux for a point source, and 
more than 90% that from a UV point source detected by 
the OM. 

Before subtracting the background flux from the 
source+background flux, we performed three count rate cal- 
ibrations, according the method described in the OM in- 
strument document The first is the deadtime correction, 
required because for a small fraction of the exposure time 
the CCD is in readout mode, and so cannot record events. 
The second calibration is for coincidence losses, which oc- 
cur when more than one photon arrives on the CCD at the 
same location and within the same frame time, so results 
in under counting. The third calibration is for the OM time 
sensitivity degradation correction. We performed these cal- 
ibrations, according to the algorithms set out in the OM 
instrument document, separately for the background and 
source+background count rates. We then subtracted the 
background count rate from the source+background count 
rate to obtain the corrected source count rate. 

Figure[T]shows the OM data points before and after cor- 
rection for aperture effects for the 17 objects. The reduced 
OM aperture does improve the alignment between the OM 
points and SDSS spectrum. This correction not only lowers 
the OM flux, but also changes the continuum shape defined 
by the OM points. Although choice of an aperture smaller 



2 URL: http://xmm2.esac.esa.int/docs/documents/CAL-TN-0019.ps.gz 

Also see the XMM-Newton User Handbook: 

http: / /xmm.esac.esa.int/external/xmm_user_support/document-ation/uhb/index.h 
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Figure 2. An example of results from SDSS spectrum fitting. The left panel shows a good fit for PG 2233+234. The black line is the 
observed spectrum, the red line is the total model spectrum. The green line represents the observed underlying continuum. The Balmer 
continuum (blue), Fell emission (light blue) and other strong emission lines (orange) are shown underneath. The right panel shows an 
example of detailed line profile fitting to the Fell subtracted region around the H/3 (upper) and Ha lines (lower) including Ha, H/3, [OIII] 
A5007/4959 doublets, [Nil] A6585/6548 doublets, Li A6708, [SII] A6717/6733 doublets, [OI] A6300/6363 doublets. In our profile fitting, 
three Gaussian components are used for H/3 and Ha, two components for [OIII] A5007, and one Gaussian for all other lines. The various 
Gaussian profiles are shown in blue, the total model is shown in red. 



than 6" will lower the OM fluxes by a larger factor, it will 
also introduce uncertainties and systematics caused by the 
PSF. Therefore we compromise by adopting a 6" diameter 
aperture. In our subsequent SED modeling we use the aper- 
ture corrected OM data. 



3 OPTICAL SPECTRAL MODELING: THE 
EMISSION LINES 

Our optical spectral modeling employs linked Ha and H/3 
profile fitting and the complete optical spectral fitting. We 
wrote the code in IDL (Interactive Data Language) v6.2, to 
perform all the optical spectral fitting. The 'MPFITEXPR' 
program from the Markwardt IDL Library is incorporated 
within our code to perform the Levenberg-Marquardt least- 
squares algorithm used to obtain the best-fit parameters. 
The SDSS spectra (stored in SDSS spSpec files) were ex- 
tracted directly from the SDSS DR7 data archive and ana- 
lyzed in IDL using our code. A detailed description of our 
spectral modeling procedures is presented in the following 
subsections. 



3.1 Profile Fitting of the Ha, H/3 and [OIIIJA5007 
Emission Lines 

Based on current AGN emission line models, there are 
thought to be stratified regions emitting different lines. 
These regions are divided somewhat arbitrarily into a 



narrow line region (NLR), a broad line region (BLR) 
and possibly an intermediate line region (ILR, e.g . 



Grup e et al. Ill999l; iHu etal1l2008l; iMei. Yuan fc Dondl2009l ; 



Zhu. Zhang fc T ang 2009). Following previous studies, we 
use several separate Gaussian profiles representing each of 
these emitting regions to model the Balmer line profiles. 

The Ha and H/3 line profiles each pose distinct difficul- 
ties for the spectral analysis. In the case of the H/3 line, the 
permitted Fell emission features (which are often strong in 
NLSls) and broad Hell 4686 line blended with the H/3 line, 
which can affect the determination of the underlying contin- 
uum and hence the H/3 line profile. For the Ha line, there is 
the problem of blending with the [Nil] A6584,6548 doublet, 
improper subtraction of which may distort Ha's intrinsic 
profile. Our approach, therefore, is to fit Ha and H/3 simul- 
taneously using the same multi Gaussian components. The 
assumed similarity between the intrinsic profiles of these two 
Balmer lines assists in deblending from other nearby emis- 
sion lines, and should yield a more robust deconvolution for 
the separate components of their profile. 

3.2 The Fell Problem 



We us e the theoretical Fell model templates of lVerner et al.l 
(2009). These include 830 energy levels and 344,035 tran- 
sitions between 2000A and 12000 A, totaling 1059 emis- 
sion lines. The predicted Fell emission depends on physical 
conditions such as microturbulence velocity and hardness 
of the radiation field, but we use the template which best 
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matc h es the observed spectrum of I Z W 1 l|Boroson fc Greer] 
1 19921 , IVeron-Cettv et alj 120041 ) i.e. the one with n H = 
10 11 cm -3 , vturb = 30 ferns -1 , Fionizing = 20.5 cm -2 s -1 . 
Detailed modelling of high signal-to-noise spectra shows 
that the Fell emission is often complex, with four ma- 
jor line systems in the case of 1 Zw 1, (one broad line 
system, two narr ow high-excitation systems and one low- 
excitation systemlVeron-Cettv et al ] |2004 IZhou et alJliobll ; 
iMei. Yuan fc Dondl2009l ~ However, for simplicity we will as- 
sume only one velocity structure and convolve this template 
with a single Lorentzian profile. 

We fit this to the actual Fell emission line features be- 
tween 5100 A and 5600 A (no other strong emission lines lie 
in this wavelength range) of the de-redshifted SDSS spectra, 
leaving the FWHM of the Lorentzian and the normalization 
of the Fell as free parameters. The resulting best-fit Fell 
model to this restricted wavelength range, was then extrap- 
olated and subtracted from the entire SDSS spectrum. A 
major benefit from subtracting the Fell features is that the 
profiles of the [OIII] A5007 lines no longer have apparent red- 
wings. This is particularly important for the NLSls, where 
the Fell emission is often strong. After subtracting Fell, we 
used either 2 or 3 Gaussian components (depending on the 
profile complexity) to fit the [OIII] A5007 line. 



3.3 Deconvolution of the Balmer Lines 

After fitting the [OIII] A5007 line, we start to fit the Ha and 
H/3 line profiles simultaneously. Following previous studies 
we consider a simplified picture in which the Balmer lines 
have three principal components, namely a narrow compo- 
nent (from the NLR), an intermediate component (from a 
transition region ILR between the NLR and BLR or fro m 
the inner edge of dusty torus (|Zhu. Zhang fc TaneH2009h l. 
and a broad component (from the BLR). The intermediate 
and broad components are both represented by a Gaussian 
profile, whereas the narrow component is assumed to be sim- 
ilar to that of [OIII] A5007. Since we do not know whether 
or not the Balmer decrements are the same in these differ- 
ent emitting zones, the relative strengths of different line 
components were not fixed, but their FWHM and relative 
velocity were both kept the same. The [OIII] A4959 line was 
set at 1/3 that of [OIII] A5007 from atomic physics. The [Nil] 
A6584,6548 line doublet were also fixed to the [OIII] A5007 
line profile. For simplicity, the [SII] A6733,6717 doublet, [OI] 
A6300,6363 doublet and Li 6708 were all fitted with a single 
Gaussian profile separately, because they are all relatively 
weak lines and do not severely blend with Balmer lines. 

In order to separate the narrow component of the 
Balmer lines from the other components as accurately as 
possible, particularly for NLSls and some broad line ob- 
jects which lack clear narrow line profiles, we applied the 
following four different fitting methods: 

1. The profile of the narrow component is held the same as 
the entire [OIII] A5007 profile; and the normalization of each 
component in the Ha and H/3 lines are left as free parame- 
ters; 

2. Only the central narrow component of the [OIII] A5007 
profile is used to define the profile of the Balmer narrow 
component, and of the [Nil] A6585,6550 doublet; the nor- 
malization of each component in the Ha and H/3 lines are 
free parameters; 



3. The shape of the narrow component is held the same 
as the entire [OIII] A5007 profile, and also the normaliza- 
tion of the H/3 line narrow component is set to be 10% of 
[OIII] A5007, t his ratio being an average for the NLR in typ- 
ical Seyfert Is (|Osterbrock fc Poggdfl985l ; iLeighlvll 19991 ); all 
other components have their normalizations as free param- 
eters; 

4. All conditions are the same as in method 3, except that 
the Balmer line narrow component and the [Nil] A6584,6548 
doublet adopt the central narrow Gaussian component of the 
[OIII] A5007 line. 

We applied each of the above fitting methods to every 
object in our sample, and then compared the results. For 
those objects with clear narrow components to their Balmer 
lines, we used the best fitting result from method 1 and 2. 
For the other objects whose narrow components were not 
clearly defined or even visible, we adopted method 3 and 
method 4, unless method 1 or 2 gave much better fitting 
results. Figure [5] right panel shows an example of our fitting. 
Results for the whole sample are shown in Figure |A"T1 

After obtaining the best-fit parameters, we used the 
intermediate and broad components to reconstruct the 
narrow-line subtracted H/3 line profile, and then measured 
the FWHM from this model. The rationale for using this 
method, instead of directly measuring the FWHM of the 
H/3 line from the data, is because for low signal/noise line 
profiles direct measurement of FWHM can lead to large un- 
certainties, whereas our profile models are not prone to lo- 
calized noise in the data. The H/3 FWHM measurements for 
each of the 51 sources, after de-convolving using the instru- 
mental resolution of 69 ferns -1 , are listed in Table [3] 



4 OPTICAL SPECTRAL MODELING 

In order to obtain the underlying continuum, we must model 
the entire SDSS spectrum so that we can remove all the 
emission lines as well as the Balmer continuum and host 
galaxy contribution. As we are now concerned with the 
broad continuum shape, we choose to refit the Fell spec- 
trum across the entire SDSS range, rather than restricting 
the fit to the Ha and H/3 line regions as discussed in the 
previous section. 

Figure [2] shows an illustrative example of our optical 
spectral fitting, and the results for each of the 51 sources 
are presented in Figure IAT1 In the following subsections we 
give further details of the components that make up these 
modeled spectra. 



4.1 Emission Lines Including Fell 

We use the models for [OIII], Ha and H/3 as derived above. 
We add to this a series of higher order Balmer lines: from 
5—5-2 (Hy) to 15— >2. We fix the line profile of these to that 
of H/3 up to 9— >2, then simply use a single Lorentzian profile 
for the rest weak higher order Balmer lines. We fix the line 
ratios for each Balmer line using the values in lOsterbrockl 
jl989h . Table 4.2, with T e between 10,000 K and 20,000 K. 
We similarly use a single Lorentizan to model the series of 
Helium lines (Hel 3187, Hel 3889, Hel 4471, Hel 5876, Hell 
3204, Hell 4686) and some other emission lines (Mgll 2798, 
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[Nelll] A3346,4326, [Oil] A3727,3730, [OI] A6302,6366, [Nil] 
A6548,6584, Li 6708, [SII] A6717,6733). 

We use the same model for the Fell emission as de- 
scribed in Section [33] However we now fit this to the entire 
SDSS wavelength range, rather than restricting the fit to 
5100-5600 A. 



4.2 The Balmer Continuum 

Another potentially significant contribut ion at shorter wave- 
length s is fr om the Balmer continu um. ICanfield fc Puetterl 
(|l98ll ) and iKwan fc Krolikl (|l98lT ) predicted the optical 
depth at the Balmer continuum edge to be less than 1, 
we use Equation [T] to model the Balmer continuum under 
the assumptions of the optically thin case and a single- 
temperature ele ctron population (also see iGrandil Il982l ; 
IWills et alJli"985h . 

F V B ° = F BB e -Hv-v BE )/(kT e) {v>Vbe) (1) 

where F BE is the flux at Balmer edge, vbe corresponds 
to the Balmer edge frequency at 3646A. T e is the electron 
temperature, h is the Planck's constant, fe is the Boltzmann's 
constant. This Balmer continuum equation is then convolved 
with a Gaussian profile to represent the real Balmer bump 
in SDSS spectra. 

There are several parameters that may slightly mod- 
ify or significantly change the shape of the Balmer contin- 
uum. It is already seen that the electron temperature T e 
appearing in Equation [T] and the optical depth can both 
change the Balmer continuum shape, but there are addi- 
tional important factors. Any intrinsic velocity dispersion 
will Doppler broaden all the Hydrogen emission features. 
Therefore a better description of the Balmer continuum can 
be obtained by convolving Equation [T] with a Gaussian pro- 
file, whose FWHM is determined by the line width of H/3 
(or other broad lines), as shown by Equation [2] where G(x) 
represents a Gaussian profile with a specific FWHM. 

/ + oo 
e -W(A*T e ) G(Ai _ A)dAi (2) 


Figure [3] shows how the Balmer continuum's shape de- 
pends on the electron temperature and velocity broaden- 
ing in Equation [2] The electron temperature modifies the 
decrease in the Balmer continuum towards shorter wave- 
lengths, but has little effect on the broadening of (Balmer 
Photo-recombination) BPR edge. On the contrary, velocity 
broadening mainly affects the shape of the BPR edge, but 
the emission longward of 3646A is still very weak compared 
to the emission blueward of the BPR edge, i.e. the BPR edge 
is still sharp. 

We initially applied Equation[2]to fit the Balmer contin- 
uum bump below 4000A in the SDSS spectra. We assumed 
the velocity profile for the convolution was a Gaussian with 
its FWHM determined from the H/3 line profile, and the 
wavelength of the position of the BPR edge was taken as 
the laboratory wavelength of 3646A. However, this model 
did not provide an acceptable fit, for example see the model 
shown by the blue line in Figure [4] It appears that the ob- 
served spectrum requires a model with either a more ex- 
tended wing redward of the BPR edge, or a BPR edge that 
shifts to longer wavelength than 3646A. However, additional 
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Figure 3. The Balmer continuum models of lGrandU d 1982T ) . The 
upper panel shows the dependence of the model on the electron 
temperature. The lower panel shows the dependence of the model 
on the FWHM of the convolved Lorentzian profile. 



velocity broadening should affect both the Balmer contin- 
uum and Balmer emission lines equally, as they are produced 
from the same material, although the multiple components 
present in the line make this difficult to constrain. 

One way the wavelength of the edge may be shifted 
without affecting the lines is via density (collisional, or 
Stark) broadening (e.g. IPigarov et al"1ll99fet ). Multiple colli- 
sions disturb the outer energy levels, leading to an effective 
n m ax for the highest bound level <C oo, i.e. lowering the 
effective ionization potential. We set the edge position and 
the FWHM as free parameters, and let the observed spectral 
shape determine their best fit values. The red line shown 
in Figure [4] represents a good fit, obtained with FWHM 
of 6000 ferns" 1 and the BPR edge wavelength of 3746 A, 
which implies n max ~ 12. The theoretical n max can be de- 
termined by the plasma de nsity N e and t emperature T e as 
n max = 2 x 10 4 (T e /iV e ) 1/4 (|Mihaladll97g| ). so for a typical 
temperature of 10 4 ~ 5 K, the required density is 7 x 10 16 ~ 17 
cm -3 . Such high density is not generally associated with the 
BLR clouds, and may give support to models where the low 
ionization BLR is from the illum inated accretion disc (e.g. 
ICollin-Souffrin fc Dumontlll990h . However, any reliable es- 
timation of the density would require more accurate sub- 
traction of other optical components such as the Fell line 
blends and many other non-hydrogen emission lines, which 
is not the focus of this paper. Nonetheless, this remains an 
interesting problem which is worthy of further study. 

Yet another issue in modeling the Balmer continuum 
is how to quantify the the total intensity of this continuum 
component, especially when there is limited spectral cov- 
erage bellow 4000A, which makes it difficult to define the 
overall shape. The theoretical flux ratio between the Balmer 
continuum and the H/3 li ne under case B c onditions can be 
expressed by Equation 131 (jWills et al.lll985h . 

I(Bac) /I(H/3) = 3.95 Tl A (3) 
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Figure 4. An expanded view of the region around the BPR edge 
in PG 1427+480. The blue and dashed lines represent the Balmcr 
continuum model superposed on the underlying disc continuum 
(green solid line) using standard parameters (blue dash), and also 
a set of best fit parameters (red dash line). The red and blue solid 
lines are models of the total optical spectrum, including the cor- 
responding Balmer continuum components and plus other com- 
ponents described in the text. The observed spectrum is shown 
in black. 



but other theoretical calculations of photonionization mod- 
els show that by varying the Balmer optical depth, elec- 
tron temperature and electron number density, this can 
resu lt in very different v alues of I(Bac)/I(H/3). For exam- 
ple, ICanfield fc Puetterl i|l98ll) 's calculation resulted in a 
I(Bac)/Ha range of 0.05~10, iKwan fc Krolikl |l98ll ) sug- 
gested I(Bac)/I(H/3) = 1.6~15, and other theoretical work 
also confirmed a large range in flux ratios (|Puetter fc Levanl 
1 19821 ; iKwanll 19841 ; iHubbard fc Puetterl [19851). The observed 
ranges in I(Bac)/I(H/J) are also large. ICanfield fc Puetterl 
(1 1981T) showed an observed range of 0.5~3 for I(Bac)/I(Ha). 
IWills et all (|l985l ) observed 9 intermediate redshift QSOs 
whose I(Bac)/I(H/3) ranges from 4.65~9.5. Thus we were 
unable to constrain the intensity of the whole Balmer contin- 
uum by using a standard flux ratio fixed to the other Balmer 
emission lines. As a result, we must rely on the shape of the 
observed Balmer bump, and then adopt the model's best fit 
parameters. 

However, this limitation in defining the Balmer bump 
introduces uncertainties in modeling the underlying contin- 
uum, because over-subtraction of the Balmer bump will de- 
press the slope of the remaining underlying continuum, and 
vice-versa. In the course of the broadband SED fitting de- 
scribed in section [5] we found that the temperature of ac- 
cretion disc (determined by black hole mass) is sensitive to 
the slope of optical continuum, unless the continuum slope 
is in the opposite sense to that of the accretion disc model 
and thus can not be fitted, or there are OM points provid- 
ing stronger constraints. We also found that a flatter optical 
continuum may lead to a lower best-fit black hole mass, al- 
though this also depends on other factors. Therefore, the 
subtraction of the Balmer continuum can have an impact 
on the modeling of broadband SED and the best-fit black 
hole mass. The influence of this depends on the relative im- 
portance of other SED restrictions. This is the reason why 



the Balmer continuum must be carefully modeled and sub- 
tracted. 



4.3 The Intrinsic Underlying Continuum 

Our basic assumption is that the residual optical spectrum, 
after subtraction of the Balmer continuum, Fell emission 
and other emission lines mentioned previously, arises mainly 
from the accretion disc emission. As a reasonable approxi- 
mation over a limited wavelength range we use a powerlaw 
of the following form to fit the underlying continuum, 

F(X) = Ci-(A/5100A)~ C2 (4) 

The powerlaw approximation for the optical underlying 
disc continuum is also widely adopted in previous and 
recent AGN optical spectral studies, (e.g. Grandi _ 19821 ; 
iTsuzuki et alll2006l : IZhou et al.ll200rj ; lLandt et all [20 111). 

We model the dust reddening using the lSeatonl (|l979l )'s 
HOOA to 10000A reddening curve, and we apply this to 
the overall model, i.e. emission lines, Balmer continuum and 
the disc continuu m. There ar e also other reddening curves 
available su ch as Fitzpatrickl (1 19861 ) for the Large Magel - 
lanic Cloud, IPrevot et all l| 19841 ) a ndlBouchet et all (Il985l ) 
for the Small Magellanic Cloud and lCalzetti et all (|2000h for 
starburst galaxies, but over the wavelength range of 2500A 
to 10000A, the d ifference betwe e n thes e reddening curves is 
small, except for ICalzetti et all {2000 )'s curve which is ap- 
propriate for starburst galaxies, and is thus not applicable 
for our AGN sample. 



4.4 The Host Galaxy Contribution 

Many previous studies on AGN's optical/infrared spectra 
have adopted a powerlaw as a reasonable approximation 
for the accretion disc c o ntinuum blueward o f lfim (e.g. 
iMei. Yuan fc Pond [20091 ; iBian fc Huanel |2010| ). but these 
studies also needed to include additional contributions from 
the host galaxy and emission from the dusty torus to account 
for the extra continuu m emission at long wavelengths of the 
optical spectrum (e.g. iKinnev et all Il996l ; iMannucci et all 
l200ll ; lLandt et alll201ll l. In our work we have also identified 
an inconsistency between the 3000A~8000A spectral shape 
and a single powerlaw shape (i.e. the flat optical spectrum 
problem discussed in Section I5.3.2[) . The blue end of the op- 
tical spectrum, presumed to arise from a standard accretion 
disc, often shows a steeper spectral slope than the red end. 

However, in our sample we found evidences suggesting 
only a weak if any, contribution from the host galaxy. For 
example, the optical spectra of our sample do not show the 
strong curvature characteristic of the presence of a stellar 
component in a host galaxy. Furthermore, the good qual- 
ity optical spectra do not exhibit stellar absorption features 
(see Section 15 . 3 . 2 1 and Figure [SJ. In fact the 3" diameter fibre 
used to obtain the SDSS spectra also helps to reduce the con- 
tribution of stellar emission from a host galaxy, particularly 
for nearby sources in our sample such as KUG 1031+398. 
These evidences argues against the possibility that the red 
optical continuum is primarily dominated by host galaxy 
emission. In fact, it is possible that the observed additional 
component arises due to emiss ion from the outer regions of 
a standard accretion disc (e.g. ISoria fc Puchnarewical2002l ; 
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ICollin fc Kawaguchill2004l ; lHao et all2010l ). The existence of 
such an additional red optical continuum component reduces 
the consistency of a powerlaw fit to the optical spectra. 

4.5 The Optical Spectrum Fitting 

Our optical spectral fitting is performed only for data blue- 
ward of 7000 A. The choice to truncate the model at 7000 A 
is made for several reasons. We wish to include Ha line in 
the spectral fitting range, and the broad wing of Ha pro- 
file sometimes extend to ~7000A (e.g. PG 1352+183, RBS 
1423, Mrk 926). There are some objects whose SDSS spec- 
tra extend only to -6700A (e.g. 2XMM J080608.0+244421, 
HS 0810+5157, 2XMM J100025. 2+015852). The choice of 
7000A, rather than a longer wavelength, is to maintain con- 
sistency of optical spectral fitting for the whole sample. 
The final reason concerns an aspect of the powerlaw fit- 
ting. We found that in some objects (e.g. PG 1115+407, 
LBQS 1228+1116, PG1352+183), a flat slope power-law un- 
der predicts the observed emission at ~7000A. Therefore, if 
we include longer wavelengths than 7000A, our powerlaw 
fitting for the standard accretion disc continuum towards 
the blue optical spectra would be biased by other contin- 
uum emission at these longer wavelengths, and so affect the 
broadband SED fitting. Consequently, we chose to truncate 
our optical spectral fitting at 7000A. 

However, we still cannot be sure that the underlying 
continuum is totally free from other non-disc continuum 
components. So after completing the fitting procedure, we 
then checked the spectral fitting status within two narrow 
wavebands, i.e. 4400A - 4800A and 5100A - 5600A. Emis- 
sion features if present in these two wavebands are mainly 
from Fell emission, and the underlying continua of these two 
wavebands should be totally dominated by the accretion disc 
emission. Assuming that the Fell emission lines within these 
two wavebands have similar relative intensity ratios as in the 
Fell template described in Section \3. 21 the best-fit underly- 
ing powerlaw plus Fell emission model should have good 
fitting status in both of these two wavebands. In general, 
the best-fit model derived from the full optical spectrum fit 
also gives reasonably good fitting status in both of these two 
narrow wavebands. However, in some cases the model over- 
predicted the flux in 5100A - 5600A but under-predicted the 
flux in 4400 A - 4800 A, so that we should slightly increase 
the slope of powerlaw to produce better spectral fitting in 
these two wavebands. We adopted these parameter values 
in preference to those directly from the full spectrum fit, 
as they should be more immune to problems such as host 
galaxy or hot dust contamination. 



5 THE BROADBAND SED MODELING 
5.1 Data Preparation 

For each object we extracted the original data files (ODFs) 
and the pipeline products (PPS) from XMM-Newton Sci- 
ence Archive (XSA) 0. In the following data reduction pro- 
cess, tasks from XMM-Newton Science Analysis System 

3 http: //xmm.esac.esa.int/external/xmm_data_acc/xsa/index. 
shtml 



(SAS) v7.1.0 were used. First, EPCHAIN/EMCHAIN tasks were 
used to extract events unless the events files had already 
been extracted for each exposure by PPS. Then ESPFILT 
task was used to define background Good Time Intervals 
(GTIs) that are free from flares. In each available EPIC im- 
age, a 45"radius circle was used to extract the source region, 
and an annulus centered on the source with inner and outer 
radii of 60"and 120"was used to define the background re- 
gion. For other sources listed in the region files of PPS that 
are included in these regions, these were subtracted using the 
default radii generated by PPS, which scaled with the source 
brightness. Then the GIT filter, source and background re- 
gion filters were applied to the corresponding events files 
to produce a set of source and background events files. We 
only accepted photons with quality flag =0 and pattern 0~4. 
The EPATPL0T task was then used to check for pile-up ef- 
fects. When pile-up was detected, an annulus with inner and 
outer radii of 12"and 45"was used instead of the previous 
45"radius circle to define the source region. Then source 
events files were reproduced using the new source region 
filter. Source and background spectra were extracted from 
these events files for each available EPIC exposure. Tasks 
RMFGEN/ARFGEN were used to produce response matrices and 
auxiliary files for the source spectra. These final spectra were 
grouped with a minimum of 25 counts per bin using the 
GRPPHA v3.0.1 tool for spectral fitting in Xspec v 11.3.2. 
To prepare the OM data, the om_filter_default.pi file and 
all response files for the V,B,U, UVW1, UVM2, UVW2 fil- 
ters were downloaded from the OM response file directory 
in HEASARC ArchivcQ. We then checked the OM source 
list file for each object to see if there were any available OM 
count rates. Each count rate and its associated error were 
entered into the om_filter_default.pi file and then combined 
with the response file of the corresponding OM filter, again 
by using the GRPPHA tool to produce OM data that could be 
used in Xspec. 

Finally, the XMM-Newton EPIC spectra are combined 
with the aperture corrected OM photometric points, and the 
optical continuum points produced from the optical under- 
lying continuum (obtained from the full optical spectrum 
fitting) using FLX2XSP tool. From these data we constructed 
a broadband nuclear SED of each AGN. There is a ubiqui- 
tous data gap in the far UV region which is due to photo- 
electric absorption by Galactic gas. Unfortunately, in most 
cases of low-redshift AGN, their intrinsic SED also peaks in 
this very UV region, and so this unobservable energy band 
often conceals a large portion of the bolometric luminosity. 
In order to account for this, and to estimate the bolomet- 
ric luminosity, we fit the X-ray and UV/optical continua all 
toget her using a new broadband SED model |Done et al.l 
l201ll . Xspec model: optxagn). We then calculate the bolo- 
metric luminosity by summing up the integrated emission 
using the best-fit parameters obtained for each continuum 
component. 

5.2 The Broadband SED Model 

A standard interpretation of the broadband SED is that 
the emission is dominated by a multi-temperature accretion 

4 http:/ /heasarc. gsfc.nasa.gov/FTP/xmm/data/rosponses/om/ 
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Table 2. Broadband SED Fitting Parameters, and Model Outputs (LboZ > fd, fa fp)- ID: object number, the same as Tablc[T] 
^H,gal and ^H,int'- the fixed galactic and free intrinsic neutral hydrogen column densities in 10 20 cm -2 ; V pow : the powerlaw 
component's slope in the SED fitting, (*) denotes the objects whose powerlaw slopes hit the uplimit of 2.2 and were fixed 
there; Fpl: the fraction of powerlaw component in the total reprocessed disc emission; R cor : corona (truncation) radius in unit 
of Gravitational radii (r 3 ) within which all disc emission is reprocessed into the Comptonisation and powerlaw components; 
T e : temperature of the Compton up-scattering electron population; Tau: optical depth of the Comptonisation component; 
log(Mg£f): the best-fit black hole mass; log(M): total mass accretion rate; L;, ;: bolometric luminosity integrated from O.OOf keV 
to 100 keV; f^, f c , f p : luminosity fractions of disc emission, soft Comptonisation and hard X-ray Compotonisation components 
in the bolometric luminosity; \ 2 '- * ne reduced \ 2 of the broadband SED fitting. 
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Table 3. Broadband SED Key Parameters. ID: object number, the same as Table [T] r2_iofceV ; the slope of the single 
powerlaw fitted to 2-10 keV spectrum. ~L2-iokeV- 2-10 keV luminosity (in 10 44 erg s _1 ); «2— lOfcev: the 2-10keV bolometric 
correction coefficient; j^: ^ e monochromatic luminosity at 2500A (in 10 43 erg s — 1 ); vYj^knV'- the monochromatic 

luminosity at 2keV (in 10 43 erg s — 1 ); a ox : the optical X-ray spectral index; AL5100: the monochromatic luminosity at 
5100A (in 10 44 erg s" 1 ); K5100: the 5100A bolometric correction coefficient; FWHMjja: the narrow component subtracted 
H/3 FWHM; L bol /L Edd : the Eddington Ratio. 
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disc c omponent which peaks in the UV (e.g.fGicrliriski ~t al.l 
Il999l . Xspec model: diskpn). This produces the seed photons 
for Compton up-scattering by a hot, optically thin elec- 
tron population within a corona situated above the disc, 
resulting in a power law component above 2 keV (e.g. 
Haardt fe Maraschil Il99ll ; IZdziarski. Poutanen fe Johnson! 
20001 . Xspec model: bknpl). However, the X-ray data 



clearly show that there is yet another component which 
rises below 1 keV in almost all high mass accretion 
rate AGNs. The ubiquity of this component can be 
seen, for exa mple, in the compilation of AGN SEDs 
presented in iMiddleton. Done fc Gierliriskil (|2007l ). and 
one of the strongest cases is th e NLS1 RE J1034+396 
i|Casebeer, Leighlv fc Baronl 120061; IMiddleton et al l 120091 ) 
and RX J0136.9-3510 ( jjin et al.1 I2009T I. The origin of 
this so-called soft X-ray excess is still unclear (e .g. 



Gierlinski fc Donell2004 ICrummv et al.ll2006l : iTurner et all 
20071 : iMiller. Turner fc Reeved^OOsf TaJid so some previous 



broadband SED modeling studies have explicitly excluded 
data below 1 keV. An obvious consequence is that in such 
studies a soft excess component cannot influence the mod- 
els, so making it possible to fit the dat a using just a disc and 
(brok en) power law continuum (VF07; IVasudevan fc Fabianl 
2009). However, in our current study we include all of the 
data, and so we require a self-consistent model which incor- 
porates this soft component. 

Whatever the true origin of the soft X-ray excess, the 
simplest model which can phenomenologically fit its shape 
is the optically thick, low temperature thermal Comptonisa- 
tion model (compTT). But the observed data are used to con- 
strain the three separate components, discpn + compTT + 
bknpl, which is generally problematic given the gap in spec- 
tral coverage between the UV and soft X-ray regions caused 
by interstellar absorption. So instead, we combine these 
three components together using a local model in XSPEC, 
assuming that they are all ultimately powered by gravita- 
tional energy released in accretion. A complete description 
of this model can be found in the Xspec website 5 and is also 
given m iDone et al.1 (|201lT ). It is in essence a faster version 
of the models recently applied to black h ole binary spec- 
tra o bserved close to their Eddington limit |Done fc Kubotal 
2006) and to the (p ossibly super Eddington) Ultra- 
Luminous X-ray sources l|Gladstone. Roberts fc Donell2009l ; 
IMiddleton fc Dondlioiol ). thus this model is more appropri- 
ate for fitting a medium sized sample of objects. A com- 
prehe nsive compa rison with t h e mod el of IDone fc Kubotal 
( 2006) is given in IDone et al.1 (|201ll ). To make this paper 
self contained we give a brief synopsis of the model. We as- 
sume that the gravitational energy released in the disc at 
each radius is emitted as a blackbody only down to a given 
radius, R CO rona- Below this radius, we further assume that 
the energy can no longer be completely thermalised, and is 
distributed between the soft excess component and the high 
energy tail. Thus the model includes all three components 
which are known to contribute to AGN SED in a self con- 
sistent way. As such it represents an improvement on the 
fits in VF07 in several respects, by including the soft excess 
and by requiring ener gy conservation, and it improves on 
IDone fc Kubotal |2006l ) by including the power law tail. 

In our SED fitting, the optical/UV data constrains the 
mass accretion rate through the outer disc, provided we have 
an estimate of the black hole mass. We constrain this by our 



analysis of the H/3 emission line profile. The main difference 
from previous studies based on non-reverberation samples is 
that we do not directly use the FWHM of the H/3 profile to 
derive the black hole mass. Rather, we use the FWHM of the 
intermediate and broad line component determined from the 
emission line fitting results pr esented in Section 13. ll These 
are then used in Equation [5] (|Woo fc Urrvl [20021 and refer- 
ences therein) to derive the black hole mass limits required 
for the SED fitting: 



Mbh = 4.8nx\ XL ^ 5W0A ) rFWHM 2 



(5) 



where La(5100A) is measured directly from the SDSS spec- 
tra. The rms difference between the black hole masses from 
this equation and from the reverberation mapping study is 
~0.5 dex. Thus we also adopted any best-fit values that fell 
below the original lower limit (which was set by FWHM of 
the intermediate component) by less than 0.5 dex. With this 
method, the best-fit black hole mass found by SED fitting 
is always consistent with the prediction from the H/3 pro- 
file. Section T6 . 5 1 discusses the differences between the best-fit 
black hole masses and those estimated using other methods. 

Once the black hole mass is constrained, the optical 
data then sets the mass accretion rate M, and hence the 
total energy available is determined by the accretion effi- 
ciency. We assume a stress-free (Novikov-Thorne) emissiv- 
ity for a Schwarzschild black hole, i.e. an overall efficiency 
of 0.057 for Ri n = 6R g . Thus the total luminosity of the 
soft excess and power law is 0.057Mc 2 (l — Rin/ Rcorona)- 
This constrains the model in the unobservable EUV region, 
with the input free parameter R CO rona setting the model 
output of the luminosity ratio between the standard disc 
emission and Comptonisation components. The upper limit 
of Rcorona is set to be 100 R g , which corresponds to 81% re- 
leased accretion disc energy. This upper limit is based on the 
requirement that the seed photons should be up-scattered 
l|Done et al.ll201ll ). We assume that both the Comptonisa- 
tion components scatter seed photons from the accretion 
disc with temperature corresponding to R CO rona- The other 
model input parameters are; the temperature (kT e ) and op- 
tical depth (r) of the soft Comptonisation component which 
are determined by the shape of the soft X-ray excess, the 
spectral index (r) of the hard X-ray Comptonisation that 
produces the 2-10 keV power law, with electron tempera- 
ture fixed at 100 keV. The model output f p i represents the 
fraction of the non-thermalised accretion energy (i.e. given 
by the luminosity originating from the region of Rcorona to 
Rin), which is emitted in the hard X-ray Comptonisation. 

We also included two sets of corrections for attenuation 
(reddening-wabs) , to account for the line of sight Galactic 
absorption and for the absorption intrinsic to each source, 
the latter is redshifted (zred and zwabs in Xspec). The 
Galactic HI column d ensity is fixed at the value taken from 
iKalberla et al.l l|2005l ). but the intrinsic HI column density 
is left as a free parameter. The standard du st to gas con - 
version formula of E(B-V) = 1.7 x 10" 22 iV H l|Besselll Il99ll ) 
is used for both Galactic and intrinsic reddening. We set 
the initial value of the powerlaw photon index to be that 
of the photon index in 2-10 keV energy band, but it can 
vary during the fitting process. However, we set an upper 
limit of 2.2 for the powerlaw photon index, not only because 
the photon index is < 2.2 for the majority of Type 1 AGNs 
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(|Middleton. Done fc Gierliriskill2007l ). but also because oth- 
erwise the much higher signal-to-noise in the soft excess in 
some observed spectra can artificially steepen the hard X- 
ray powerlaw and result in unphyiscal best-fit models. 

All free parameters used in the broadband SED fitting 
are listed in Table [2] For completeness, we also explicitly 
calculate the fraction of the total luminosity carried by each 
component of the model (i.e. disc: /<j; soft Comptonization: 
f c ; hard X-ray Comptonization: f p ) from the model fit pa- 
rameters Rcor and F p ; (see Table [2| 0. Table [3] lists the 
important characteristic parameters. The main uncertainty 
in these parameters, especially the black hole mass, is dom- 
inated by other systematic uncertainties introduced by the 
observational data, model assumptions (e.g. the assumption 
of a non-spinning black hole and the inclination dependence 
of the disc emission) and the analysis methods involved. 
Therefore the parameter fitting uncertainties which are often 
less than 10%, are not significant in comparison, and thus 
are not listed. The statistical properties of these parameters 
are discussed in section [6] 



5.3 Problems in The SED Fitting 

We further discuss two problems we encountered during the 
fitting procedure in the following subsections. The first prob- 
lem is the discrepancy between the OM and SDSS contin- 
uum points (mentioned in Section [23]). The second problem 
is that of the observed flat optical continuum, whose shape 
cannot be accounted for in our SED model (mentioned in 
Section |4~4)) . 



5. 3. 1 The discrepancy between the OM photometry and 
the SDSS continuum 

There remains a significant discrepancy between many of 
the OM and SDSS continuum points, even after applying the 
aperture correction discussed in Section [272] (see Figure [AT]) , 
The OM points often appear above (brighter) the extrapo- 
lation of the SDSS continuum to the OM wavelengths. We 
identify three possible reasons for this discrepancy: 

(1) Remaining aperture effects: There is an aperture 
difference between the SDSS fibres (3"diameter) and the 
OM apertures we used (6"diameter). Clearly the OM points 
will still include more host galaxy starlight than the SDSS 
points, and so will appear above the SDSS spectrum. 

(2) Contamination from emission lines: The wavelength 
ranges for each OM filter (over which the effective trans- 
mission is greater than 10% of the peak effective trans- 
mission) are as follows: UVW2 1805-2454A, UVM2 1970- 
2675A, UVW1 2410-3565A, U 3030-3890A, B 3815-4910A, 
V 5020-5870A. We exclude the contribution from strong op- 
tical emission lines within the OM U, B, V bandpass (and 
also the Balmer continuum contribution in U band) by us- 
ing the best-fit optical underlying continuum which excludes 
such features from the SDSS spectral fitting. In fact, this 



5 a full description of the model parameters can be found on the 
Xspec web page: 

http://hcasarc.nasa.gov/xanadu/xspec/models/optxagn.html 



was an important initial motivation of the study, i.e. to ob- 
tain more accurate estimates of the true underlying contin- 
uum rather than simply to use the SDSS 'ugriz' photometric 
data. Inclusion of strong emission lines within these photo- 
metric data would result in over-estimation of the optical 
continuum, and so compromise our aim to study the shape 
of the optical underlying continuum. This is an important 
spectral characteristic used to constrain the accretion disc 
component in the SED fitting (see also the discussion in 
section 5.1.2). There are some strong emission lines within 
the UV bandpasses such as Lya, CIV 1549, CIII 1909 and 
Mgll 2798, whose fluxes are not available from SDSS spec- 
trum. Accurate subtraction of these line fluxes for each ob- 
ject would require new UV spectroscopy. We conclude that 
inclusion of emission line flux within the OM photometric 
points may account for some of the observed discrepancy. 

(3) Intrinsic source variability: AGN are well known to 
be variable across their SEDs. In general there is a signif- 
icant time difference between acquistion of the SDSS and 
OM-UV data, so intrinsic variation may contribute to any 
observed discrepancy. Mrk 110 is the most extreme exam- 
ple of this phenomena in our sample, as its SDSS spectrum 
has a very large dis crepancy compared with the OM data. 
The recent paper bv lLandt etail (|201ll ) gives another set of 
optical spectra for Mrk 110, which is more consistent with 
our best-fit model. It shows that the inclusion of OM data 
is useful to help identify cases such as this. As an additional 
test for variability, we assembled all available GALEX data 
for our sample. We find that 43 objects in our sample have 
GALEX data. Using a GALEX aperture of 12", which is 
limited by the PSF and which is also similar to the UV OM 
apertures, we compare these values with the SED model. 
The ratio of the GALEX data and our SED model within 
the same bandpass differ by less than a factor of 2 for the 
majority of our sample, and significantly the flux ratio dis- 
tribution is almost symmetric and is centered close to unity. 
This suggests that the non-simultaneous OM and SDSS data 
is not likely to be a major impediment to our modeling. 

In effect, these three factors will merge together to pro- 
duce the observed discrepancy between the SDSS and OM 
data. Since the combined effects of Point (1) and (2) which 
will add flux and generally be greater than that caused by 
optical/UV variability as sh own by previous long term rever- 
berat ion mapping studies (|Giveon et al ] |l999l :l Kaspi et al.l 
2000J), we should treat the OM points included in our SED 
modeling as upper limits when interpreting the results of 
our modeling. Indeed, the 90% confidence uncertainties in 
the BH masses derived directly from the Xspec fitting are 
almost certainly small compared with the systematic errors 
introduced by the above uncertainties. 

5.3.2 The observed flat optical continuum 

A related problem in our fitting is about the SDSS contin- 
uum shape. For some AGNs, their SDSS continuum data 
points exhibit a very different spectral slope from that of 
the SED model. This cannot be reconciled by adjusting the 
parameters of the accretion disc model, and thus implies 
the presence of an additional component at longer optical 
wavelengths, which flattens compared with that predicted 
by the accretion disc models. One obvious explanation for 
this flux excess is the contribution from the host galaxy. In 
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Figure 5. A comparison between the results of two subtractions of host galaxy contribution. 2XMM J112328+052823 (Fig-al and Fig-a2) 
shows an underlying continuum that more closely resembles a disc continuum (solid green line in Fig-al) after modelling and subtracting 
the host galaxy contribution (light blue spectrum in Fig-al). The left panel of Fig-a2 shows the original broadband SED fitting without 
subtracting the host galaxy contribution. The dashed green line shows the modelled accretion disc emission in the best-fit SED. The 
inserted panel shows a magnification of the fit in the optical/UV region, where a big discrepancy exists between the SDSS data and 
best-fit SED model. The right panel of Fig-a2 is the new SED fit using the new underlying disc continuum (shown as solid green line in 
Fig-al) after subtracting the host galaxy contribution. The new fit is improved in the optical region compared with the previous results 
in the left panel of Fig-a2. In contrast to the above example, PG 1415+451 (Fig-bl and Fig-b2) has little host galaxy contribution in 
the SDSS optical spectrum (see the light blue component in Fig-bl), and its broadband SED fitting in the optical region remains poor 
re gardless of the amou nt of host galaxy subtraction applied (see the two panels in Fig-b2). The spectral template for Elliptical galaxies 
in lKinnev et al.l jl996) was used in both cases since their host galaxies both have elliptical morphologies in SDSS image. 



late type host galaxies such as elliptical and SO galaxies, 
emission from their old stellar population s peaks at near 
infrared wavelengths. iKinnev et al.l (|l996T ) combined spec- 
tra of quiescent galaxies and constructed an average spec- 
tral template for each morphological type, including bulge, 
elliptical, SO, Sa, Sb, Sc and starburst galaxies. For some 
objects in our sample with high S/N SDSS spectra which 
show at least marginal stellar absorption features, we have 
added the corre sponding type of hos t galaxy spectral tem- 
plate taken from lKinnev et all (1 19961 ). into the overall SDSS 
spectral fitting. This revised underlying continuum in the 
optical, and was then used in the broadband SED fitting. 
We are then able to compare it with the original fit, to see 
how the subtraction of a stellar population template effects 
the overall SED fitting. 

Figure [5] shows two examples. The first is 2XMM 



J112328.0+052823, in which after subtracting the host 
galaxy component, the observed optical continuum is closer 
to the slope of the SED model. However, the results for 
PG1415+451 in Figure [5] lower panel imply that its host 
galaxy cannot be the origin of the flat optical spectrum. 
The reason is that its optical spectrum does not show any 
strong stellar absorption features. This means that the maxi- 
mum amount of host galaxy contribution is small, and so and 
there remains a substantial inconsistency in the slope versus 
the SED model. In addition to 2XMM J112328.0+052823 
above, only Mrkl018 and 2XMM J125553.0+272405 show 
clear stellar absorption features. Also the 3"diameter fibre 
excludes much of the host galaxy component at these red- 
shifts. Therefore, on these general grounds we conclude that 
host galaxy contamination is small for most sources in our 
sample, and consequently cannot fully account for the ob- 
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served flat optical continuum. Additional support for this 
view comes from good correlations between the X-ray com- 
ponents and the red optical continuum, suggesting that this 
extr a optical flux is likely related to the intrinsic activity 
(e.g. Soria fc Puchnarewiczll2002l: I Collin fe Kawaguchill200l 



lHao et al.ll2010l : lLandt et al]l201ll ) 



6 STATISTICAL PROPERTIES OF THE 
SAMPLE 

Histograms of data on our sample are shown in Figure [6l 
Figure [8] and Figure [9] including redshift, HI column density, 
optical and X-ray modeling parameters etc. The red region 
in the histograms show the distributions for the 12 NLSls 
in our sample. It is clear that NLSls are distinct among the 
whole sample in several respects. 



6.1 General Properties 

Figure [6] shows some basic properties of our sample which 
are not model dependent: 

(1) . Redshift: the sample's redshift ranges from 0.031 (Mrk 
493) to 0.377 (HS 0810+5157). The NLSls are found mainly 
at lower redshifts, with < z > n = 0.12 compared to the 
< z >„~ 0.19 for the BLSls. For comparison we see that 
the sample of VF07 has a similar redshift range, but it has 
a lower average redshift of 0.10. 

(2) . The Galactic nH: the average Galactic nH is 2.25xl0 20 . 

(3) . The photon indexes obtained from simple power law 
fits to the restricted energy range of 2-10 keV. The NLSls 
cluster on the higher photon index side, with an aver- 
age of 2.21±0.20, which differs from the sample average of 
1.92±0.25 and the BLSls' average of 1.83±0.18. This means 
that NLSls tend to have softer X-ray spectra, which is fur- 
ther confirmed in the following section on the mean SEDs. 

(4) . The X-ray continuum and 2-10 keV luminosity: this dis- 
tribution shows that NLSls have lower 2-10 keV luminosities 
in spite of their steeper slopes. We note that the VF07 sam- 
ple has a similar distribution, except for their inclusion of 
three extremely low X-ray luminosity AGN (i.e. NGC4395, 
NGC3227 and NGC6814), these objects were not included 
in our sample due to our selection criteria and/or a lack of 
SDSS spectra. 

(5) . The optical continuum luminosity at 5100 A. On aver- 
age the NLS1 have lower optical luminosities than BLS1. 
(6-8). The [OIII] A5007, Ha and R/3 emission line luminosi- 
ties. Again the NLSls have on average lower luminosities 
than BLSls. 

(9). The Balmer decrement. The average value for the whole 
sample is 3.14±0.62, and for NLSls is 3.05±0.38. This dif- 
ference is not statistically significant, but we return to the 
issue in our next paper (paper II), where we consider the 
separate components as well as the overall profile. 

6.2 Results from The Broadband SED Modeling 

Figure [7] shows properties derived from the SED fits: 
(1). The bolometric luminosity: the distribution range is be- 
tween 1.8xl0 44 erps s" 1 (Mrk 464) and 1.4xl0 47 ergs s" 1 
(PG 2233+134). There is no clear difference in the distribu- 
tion of the complete sample and the sub-set of NLSls. The 



average luminosity is Log(L(, o ;)=45.49±0.55, which is con- 
sistent with the value of 45.19+1.01 found in VF07 sample, 
except for the three extremely nearby and low luminosity 
AGNs in VF07. 

(2) . The black hole mass: using the best-fit black hole 
masses, the whole sample peaks between 1O 7 M0 and 
10 8 Mq. Equation[5]suggests that the black hole mass should 
depend on both H/3 FWHM and L5100 , and the results from 
our SED fitting suggest that NLSls with smaller Balmer 
line FWHM do indeed harbour lower mass black holes. KUG 
1034+396 has the lowest black hole mass in our sample. The 
value of 1.7x 1O 6 M0 is consistent with the estimate based on 
the first firmly detected A GN QPO (quasi perio dic oscilla- 
tion) found in this source l|Gierlihski et al.ll2008f ). A gain we 
can compare our results with those of VF07 sample. We find 
that their average black hole mass is 7.89+0.82, calculated 
using the M(Lsioo, FWHMjj^) relation. Adopting this same 
method for our sample, we find a very similar average of 
7.99+0.93. Our best-fit masses have a slightly lower average 
value of 7.83+0.64 (also see Section \6. 5 1 for a comparison of 
different estimates of black hole masses). 

(3) . The Eddington ratio: the average values are 3.21+3.07 
for NLS1 which display a wide dispersion, and 0.57+0.50 
for BLS1 and 0.93+0.85 for the whole sample. Of the eight 
objects whose Eddington ratios are above 1, six are NLS1 
galaxies, and the highest value is 14.2 (PG 2233+134). 
Clearly, NLSls tend to have larger Eddington ratios. Our 
Eddington ratio distribution is also similar to that found in 
the sample of VF07 whose average value is 0.47+0.44, ex- 
cept that their distribution has a more pronounced peak at 
~0.1. 

(4) . The a x index, is defined bet ween restframe co ntinuum 
points at 2500 A and 2 keV (see iLusso et al.ll20ld and ref- 
erences therein). The distribution for NLS1 is peaked at 
marginally higher values than for BLS1. 

(5) . The K2-10 bolometric correction, is defined as 
Lboi/L2-io (see VF07 and references therein). We find that 
NLSls have a significantly higher fraction of their bolometric 
luminosity emitted as hard X-rays than the BLSls. Com- 
pared with the VF07 sample, both distributions peak at 
K2-io=10~30, but our sample shows a smoother distribu- 
tion decreasing as K2-10 increases after ~30, and so results 
in a slightly higher average value of «2-io- 

(6) . The intrinsic nH: This distribution shows that the in- 
trinsic equivalent neutral hydrogen column densities are low 
for our sample, which is a natural consequence of our initial 
sample selection criteria. The NLSls have slightly higher in- 
trinsic absorption than BLSls, which may imply a slightly 
higher dust reddening. However the distribution of Balmer 
decrements shows no significant difference between these two 
types of AGNs. 

(7) . The temperature of the Comptonisation component 
used to describe the soft X-ray excess. This is close to 
0.2 keV in all objects, confirming the trend seen in previ- 
ous studies fo r this component to exhibit a narrow range o f 
peak energy l|Czernv et al J 120031 : iGierlihski fc Done] |2004| ) . 
The distribution peak at this energy is more marked for the 
BLS1 than for NLS1, although the small number statistics 
means that this difference cannot be considered as definitive 
for our sample. 

(8) . The optical depth of the soft excess Comptonised com- 
ponent. It is clear that this component is always optically 
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Figure 6. Distributions of our sample for different properties. In each panel the blue areas show the distribution for the whole sample, 
while the red areas show the distribution for the 12 NLSls in our sample. We note that the Ha, H/3 and [OIII] A5007 luminosities are 
based on results of line profile fitting, after subtracting the blends from other nearby emission lines (see Section |3. It . For comparison we 
also indicate the Balmer decrement value of 2.86, found under case B recombination, as shown by dashed line in the same panel. 



thick, with most objects having r~10 — 30. There is no sig- 
nificant difference in temperature or optical depth between 
the broad and narrow line objects. 

(9). There is a difference in the coronal radii distribution 
between the BLSls and NLSls. Corona radius controls the 
relative amount of power emerging from the accretion disc 
and the soft X-ray excess/hard tail. There are two peaks 
in the distribution for the broad line objects, one between 
10 and 20 R 9 (where R g = GM/c 2 ), and the other at 100 
R 9 (which is set as the upper limit of this parameter in our 
broadband SED model). By contrast these radii in NLS1 
are consistent with just the first peak. At first sight this is 
surprising, since NLS1 are expected to be those with the 
strongest soft X-ray excess. However, their similar soft ex- 
cess temperatures around 0.2 keV suggests that atomic pro- 
cesses may be significant (reflection and/or absorption from 
partially ionized material), and this may influence our fits. 
The average coronal radii are 32±26 R 9 for NLS1, 59±37 R 9 
for BLS1 and 53±36 R 9 for the whole sample. This supports 
the conclusion of VF07 that high Eddington ratio AGN have 



lower coronal fractions compared to those with low Edding- 
ton ratios. 

6.3 Balmer Line Parameter Distribution 

Figure [8] shows further details of the modeled profiles of Ha 
(first row), and H/3 (second row). 

(1). The FWHM of the broad emission profile. This is cal- 
culated from co-adding the two best fit Gaussian profiles for 
the broad and intermediate line components, and then using 
the resultant profile to determine the FWHM. This is equiv- 
alent to subtracting the narrow line core from the observed 
profile and measuring the resultant FWHM. 
Note, the NLSls by definition have Ha < 2000 km s" 1 . 
(2-3). The equivalent widths and line luminosities are again 
measured using the total broad emission line profile as 
above. The NLSls have both lower equivalent widths and 
line luminosities. 

(4). By contrast, there is no pronounced difference between 
NLSls and BLSls in their Balmer narrow line component. 
This suggests that the narrow line region is less influenced 
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Figure 7. The distribution of model dependent parameters using the same colour coding as in Figure[6] Comments on each distribution 
are given in Section 16.21 



by whatever difference in properties is responsible for the 
defining difference between NLSls and BLSls in the broad 
line region. 



6.4 The Bolometric Luminosities 

The fraction of the total luminosity contained in each com- 
ponent of the SED model is shown in Figure [9] The upper 
left panels show these fractions as a function of the bolomet- 
ric luminosity. It seems that as the bolometric luminosity in- 
creases, the disc component slightly increases in importance. 
However, the total numbers of objects at high luminosities 
is small, as seen in the upper right panels, where the frac- 
tion is multiplied by the number of objects in the bin, so we 
should be cautious about this finding. 

The lower panel shows this fraction for each of the ob- 
jects ranking from the smallest to biggest H/3 FWHM. Thus 
low rank objects have the narrowest H/3 (and hence are 
by definition NLSls). These also have the lowest black hole 
masses and highest Eddington ratios. They are more likely 
to have a smaller fraction of their total luminosity emitted 
in the soft X-ray excess component, than the BLSls. This 



Table 4. The average black hole masses, as shown in Figure [Tol 





NLS1 


BLS1 


ALL 


<M BH ,IC > 


6.58±0.49 


8.09±0.56 


7.73±0.84 


<M BHtBC > 


7.72±0.49 


9.05±0.55 


8.74±0.78 


<^BH,IC+BC > 


6.75±0.49 


8.37±0.65 


7.99±0.93 


<M BHitT > 


6.57±0.46 


7.89±0.47 


7.58±0.73 


<Mbh,fit > 


7.11±0.54 


8.05±0.48 


7.83±0.64 


<M BH ,RP > 


7.42±0.39 


8.44±0.53 


8.20±0.66 



relates to the issue of the coronal radii, see Point (9) of 
Section [6.21 There are also some BLSls which have an ap- 
parently high fraction of power in their soft X-ray excesses, 
but they may also have alternative spectral fits including 
reflection and/or absorption. 

We note that in all these plots the lower limit to the 
disc fraction of 0.19 results from setting an upper limit of 
100 R 9 for the coronal radius parameter, as mentioned in 
Section [5^21 
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Figure 8. The Balmer line parameter distributions. The first row is for Ha and the second is for H/3. We combine the intermediate and 
broad components in each Balmer line profile to form the total broad line properties, giving values of the FWHM, EW and luminosity. 
The final panel shows the luminosity distribution of the narrow component for comparison. The distributions for the 12 NLSls are 
indicated by the red regions, as in Figure [6] 
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Figure 9. The bolometric luminosity distribution for the differ- 
ent continuum components of the SED, i.e. accretion disc (green), 
Comptonisation (orange) and hard X-ray Comptonisation (blue). 
The upper left panel shows the percentage within each luminos- 
ity bin for each of these three SED components. The Upper right 
panel shows the luminosity distribution of the whole sample, with 
each bin divided into three regions according to the fractional con- 
tribution from the different components in that luminosity bin. 
The lower panel shows how the contribution from each compo- 
nent changes as a function of rank order in H/3 FWHM, after the 
narrow line component has been removed. 



6.5 The Black Hole Mass 

The black hole mass is one of the key parameters used in our 
SED fitting, and it largely determines the continuum shape 
in the optical/UV region. The masses derived from rever- 
beration mapping are considered to be the most accurate, 
but the total number of objects which have b een studied us- 
ing this technique is still relatively small (e.g . iPeterson et al.l 
|2004 iDennev et al.|[201ol ; iBentz et alJl201Ch . In the absence 
of reverberation mapping, the empirical relation between 
Mbh and H/3 linewidth and L5100 is often used as a proxy 
to estimate the black hole mass (Pe terson et al 1 120041) . A 
serious limitation of this method is that it is still not clear 
which specific measure of the H/3 profile provides the clos- 
est association with the velocity dispersion of the gas in the 
broad line region. 

There are various alternative measures of the velocity 
width used for determining the black hole mass, including 
the FWHMs of the interm ediate component (IC) an d the 
broad component (BC) (e.g. lZhu. Zhang fc Tangll2009t ) . One 
coul d also use the model independent seco nd momentum 
(e.g. IPeterson etHl |2004| ; iBian et al.l 120081 ). or more sim- 
ply the FWHM of the H/3 line after subtrac ting the narrow 
component (NC) (e.g. IPeterson et al.|[200i '). The NC sub- 
tracted FWHM and the second momentum estimates often 
lie within the range of values covered by the IC and BC 
FWHMs, except for some peculiar objects such as those with 
broad double-peaked profiles, for example UM 269. Given all 
these uncertainties we decided to adopt the the best- fit black 
hole mass obtained from the SED model, rather than sim- 
ply fixing it at a value determined from a specific linewidth 
measurement. Moreover, it is now suggested that radiation 
pressure may be important in modifying the black hole mass 
derived using the relation between Mbh and L5100 and H/3 
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Figure 10. A comparison of various methods used to derive black 
hole mass. The total distributions are shown with the 12 NLSls 
show by the red regions. The purple dashed line indicate the av- 
erage black hole mass for the whole sample. The orange and cyan 
dotted lines indicate the average masses of NLSls and BLSls, 
respectively. The average values are listed in Table [4] Values for 
individual objects are listed in Table ICll 



FWHM, especially for o bjects with high Edd ington ratios 
such as most NLSls (e.g. iMarconi et al. 

In order to compare our results with those from other 
studies, we have made various estimates of black hole masses 
for every source in our sample as follows: 

(1) M BH ,ic, M B h,bc and M bh ,ic+bc are derived using 
Equation [5] with different H/3 FWHMs obtained from our 
Balmer line fitting procedure. 

(2) M bh ,a is the black hole mass calculated f rom the second 
mome ntum of the total H/3 line profile (see iPeterson et alj 
(|2004 ) for details of the definition of 'second momentum'), 
by using Rblr ocL^oq 8 and a geometry factor of / = 
3.85. These assumptions are considered to be appropriate 
when using second m omentum as a measure of the veloc- 
ity dispersion in BLR l|Bentz et al.ll2006l ; ICollin et aLll200fj ; 
iBian et al.ll2008l ). 

(3) Msff,BP is the black hole mass corrected f or rad iation 
pressure, using equation (9) in lMarcon[ et al. (2008) with 
/ = 3.1, log(g) = 7.6. 



We compare the black hole mass distributions obtained from 
these different methods in Figure 1101 The mean values are 
listed in Table g] 

The M B h,ic and Mbh.bc represent the two extreme 
estimates of black hole masses. The Mbh,ic could still be 
influenced by contamination from a NLR component, es- 
pecially for NLSls where deconvolution of the narrow and 
broad components is very difficult. If there is a residual nar- 
row line component, it will introduce a bias that under- 
estimates black hole masses. Conversely, using Mbh.bc is 
more likely to bias towards higher black hole masses, due to 
the presence of low contrast very broad wings often seen 
in U/3 profiles. We found FWHMjc+Bc/ffH/j=1.30±0.39 
for our sample , whic h is consistent with 1.33T0.36 found 
by IBian etahl (|2008l ). This leads to slightly lower values 
of Mbh, a than Mbh,ic+bc , but these two methods both 
give black hole masses between Mbh,ic and Mbh,bc, with 
Mbh, ic+bc spanning a broader mass range. 

Our best-fit SED black hole masses (Mbh,fit) are also 
distributed between Mbh,ic and Mbh,bc, with similar av- 
erage masses as Mbh,ic+bc (a comparison is shown in Fig- 
ure [11] Panel-A) . Note that Mbh, fit is a free parameter 
in the SED fitting unless it hits the lower or upper lim- 
its set by Mbh,ic and Mbh,bc, which occasionally hap- 
pened (see Tab el [Cl) . It is clearly shown in Figure [TT] that 
the black hole masses from the SED fitting are not consis- 
tent with estimates based on either extremely narrow or ex- 
tremely broad lines. So for NLSls, the mean M B h,fit is 0.36 
dex higher than Mbh,ic+bc\ while for BLSls, the mean 
M B h,fit is 0.22 dex lower than Mbh, ic+bc- Interestingly, 
this also implies that the Mbb,f/t of NLSls may have less 
deviation from the established M-Sigma relation than that 
using the M(L5ioo, F WHMfrg) relation as shown in several 
previous studies (e.g. IWang fc Lull200ll ; IBian fc Zhaoll2004l ; 
IZhou et al.ll2006T ). 

The situ at ion may be further complicated as 
IMarconi et~aT] (|2008l ) showed that NLSls could be 
consistent with the M-cr* relation if a correction for ra- 
diation pressure is applied to black hole masses derived 
from M(L5ioo, FWHMh^). In our sample, correction for 
radiation pressure adds to the average Mbh, ic+bc by 
0.67 dex for NLS1, 0.07 dex for BLS1 and 0.21 dex for 
the whole sample. We also found a very similar mass 
distribution between Mbh,rp and Mbh, fit, except for 
an average of 0.36 dex higher in Mbh,rp- The differences 
between the average mass of NLSls and BLSls are 0.78 
dex and 0.72 dex in Mbh,rp and Mbh, fit, separately (see 
Figure [TTJ Panel-B). Therefore, if Mbh,rp can provide a 
good match to the M- Sigma relation even do wn to low mass 
NLSls as proposed by IMarconi et al.1 (|2008T ), then our SED 
determined Mbh, fit may also give similar results. This 
implies that the suggested deviation from the M-a* relation 
for NLSls may not be an intrinsic property, but rather a 
consequence of using black hole estimates based on M(Lsioo, 
FWHMjf(3) relation, which may not be approp riate for 
NLSls fe.g. iGrupe fc Mathur|[200i ; Ik omoss a 2008). 

6.6 The Average Spectral Energy Distributions 

lElvis et all (1 19941 ) constructed SED templates for both 
radio-loud and radio-quiet AGN, based on a sample of 47 
quasars between redshift 0.025 and 0.94. VF07 modeled 
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Figure 11. Correlations of best-fit black hole mass ('Mg^r-Fitting' or 'Mbjj^jj'') vs. H/3 FWHM determined black hole mass ('Mbh- 
H/3 FWHM' or 'M-bh ic+bc') an d vs. radiation pressure corrected black hole mass ('M^j^-Radiation Pressure' or 'Mgjj j;p'). Red 
points represent the 12 NLSls. The inserted panel in panel-A shows the distribution of the mass difference between M-bh,ic+bc an d 
yi-BH ,FIT > while the inserted panel in panel-B shows the distribution of the mass difference between Mbh,rp an d ^Abh,fit- Red regions 
highlight the distribution of NLSls. 



optical-to- X-ray SED for a sample of 54 AGNs with red- 
shifts between 0.001 and 0.371, and showed that the SED 
was related to Eddington ratio. They also suggested that 
K2-ioke V is well correlated with E ddington ratio. In a later 
study (|Vasudevan fc Fabian! I2009TI based on SED model- 
ing of 29 local AGNs from iPeterson et all (|2004T ), the SED 
dependence on Ed dington ratio was reinforced. Recently, 
iLusso et all (|2010h studied 545 X-ray selected type 1 AGN 
over the redshift range of 0.04 to 4.25. They computed SEDs 
at different redshifts, and investigated a ox correlations with 
other parameters such as redshift, K2-iofeeV, ^Edd etc. 

We present a mean SED for our sample which is sub- 
divided according to their H/3 FWHM. This gave three sub- 
samples, those with the narrowest lines, those with mod- 
erately broad lines, and those with very broad lines. All 
objects were de-redshifted to their local frame. First, each 
of the best-fit SEDs was divided into 450 energy bins be- 
tween 1 eV and 100 keV. For each energy bin we calculated 
the monochromatic luminosity for the sub-sample with 12 
NLSls, using their individual SED models. Then an aver- 
age value and standard deviation in each energy bin were 
calculated in logarithm space. Thus a mean SED for the 
12 NLSls was constructed. Using the same method for the 
12 moderate and 12 broadest line objects, their mean SEDs 
were produced. The total SED energy range is 1 eV to 100 
keV, but we note that only spectral ranges from 1.5-6 eV and 
0.3-10 keV are actually covered by the observational data, 
and all other ranges are based on model extrapolations. 

Obviously, limitations of our mean SEDs include the 
relatively small sample sizes composing the SEDs, and the 
redshift restriction z < 0.4. On the other hand, we have 
assembled high quality data sets of optical, UV and X-ray 
observations. The exclusion of objects with high intrinsic 



absorption in the optical/UV helps to simplify the modeling 
assumptions. Our exclusion of warm-absorber objects may 
have introduced unknown selection effects, but again this 
simplified the SED modeling. Our model of the accretion 
flow also includes more detailed physical assumptions on 
the optical-to-X-ray spectrum than in previous broadband 
SED studies. These advantages make our broadband SED 
fitting more physically plausible. Thus our mean SEDs too 
should be more reliable, especially in the unobservable far 
UV region, where often the peak of the energy is emitted. 

Figure [T3] shows the mean SEDs for the three subsets of 
our sample. We caution that there is still subst antial spec- 
tral d iversity within each subsample, and echo lElvis et al.l 
l|l994h 's warning that if AGN SEDs are simply averaged 
without considering their detailed intrinsic properties, then 
the dispersion in the resultant mean SED will be large, so 
the mean SED may lose some useful information about AGN 
properties. Nevertheless, there appears to be a clear SED 
connection with H/? FWHM. As the line width increases, 
so the big blue bump (BBB) in the UV region becomes 
weaker relative to the hard X-rays, and its peak shifts to- 
wards lower energy. Also the spectral slope at high energies 
becomes harder. 

This evolu tion in spectral shape i s simi lar to that found 
by VF07 and IVasudevan fc Fabianl (|2009T l. in which two 
mean SEDs of different mean Eddington ratio were com- 
pared. This relation might be expected since the FWHM 
and Eddington ratio are also strongly (anti) correlated in our 
sample. VF07 interpreted the spectral diversity as a scaled 
up version of the different accretion states of Galactic black 
hole binaries. The low Eddington ratio AGN could be analo- 
gous to the low/hard state in black hole binaries in having a 
weak disc giving a strong high energy tail, and the high Ed- 
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dington ratio sources are analogous to the high/soft state, in 
which the disc emission dominates. Our SED templates do 
not extend down to such low Eddington ratios as in VF07, 
but we still see a similar behaviour. 



7 SUMMARY AND CONCLUSIONS 

In this paper we presented a spectral study of 51 unob- 
scured Type 1 AGNs, including 12 NLSls. We assembled 
X-ray data from the EPIC monitor on board the XMM- 
Newton satellite, and optical data from the SDSS DR7. In 
addition we added optical/UV data from the XMM-Newton 
OM monitor when available. Our results confirm some pre- 
viously known correlations. For example, NLSls often have 
softer powerlaw fits from 2-10 keV, and have lower 2-10 keV 
luminosities. Their Ha, H/3 and [OIII]A5007 lines are also 
less luminous on average than found in BLSls. 

We use detailed models to fit the Ha and Hf3 line pro- 
files, with multi-components to deblend the narrow, inter- 
mediate and broad components by means of simultaneous 
modeling of the Fell continuum and other blended lines. We 
then use results from the H/3 line fitting to constrain the 
black hole mass. The FWHM of the intermediate and broad 
components give a lower and upper limit for the mass, re- 
spectively. This supports previous studies which find that 
NLSls tend to have lower black hole masses and higher Ed- 
dington ratios, although their bolometric luminosities are 
not significantly different from those of BLSls. 

We include the Balmer continuum and permitted iron 
emission, and extend the modeling across the entire SDSS 
spectrum in order to isolate the intrinsic optical underlying 
continuum. However, this pure optical continuum is often 
(in 32/51 objects) flatter than is predicted by the standard 
accretion disc model. This could indicate some contamina- 
tion from the host galaxy, but the lack of stellar absorption 
features in most of the SDSS spectra suggests that this can- 
not be a general explanation. Instead it seems more likely 
that there is an additional component in the optical region 
related to the AGN, which is as yet not well understood. 

We also show that the Balmer continuum is not well 
modeled if the edge wavelength is fixed at its laboratory 
value of 3646A. It is shifted redwards, and smoothed by 
more than predicted by the FWHM of the Balmer emis- 
sion lines. These effects could both be produced by density 
broadening. Potentially more detailed models of the opti- 
cal emission could employ this as a new diagnostic tool for 
studying the physical conditions e.g. electron density and 
temperature, in the innermost Balmer emitting regions. 

The optical, UV and X-ray data were fitted using a 
new broadband SED model, which assumes that the gravi- 
tational potential energy is emitted as optically thick black- 
body emission at each radius down to some specific coronal 
radius. Below this radius the remaining energy down to the 
last stable orbit is divided between a soft X-ray excess com- 
ponent and a hard X-ray tail. This energetically constrains 
the model fits in the unobservable EUV region. We construct 
the resulting SEDs for each of the sources. 

A multi-component decomposition of the broadband 
SED shows that relative contributions to the bolometric 
luminosity from the accretion disc, Compotonisation and 
powerlaw components vary among sources with different lu- 



minosity and H/3 linewidth. We find a slight increase in con- 
tribution from the accretion disc as the luminosity increases, 
but a larger sample with more sources at both low and high 
luminosities is needed to confirm this. 

Our study also supports the distinctiveness of the 
NLSls among the whole sample. We find that NLSls tend 
to have a softer 2-10 keV spectrum, lower 2-10 keV lumi- 
nosity, lower black hole mass, higher Eddington ratio and 
higher a ox index. However NLSls do not stand out from 
the whole sample in terms of their bolometric luminosity 
distribution. We estimate the corona radii for every AGN in 
our sample from the SED fitting. This shows that on aver- 
age NLSls have smaller corona radii, and correspondingly a 
smaller coronal component contribution. 

We compare the best-fit black hole masses with those 
corrected for radiation pressure, and other estimates of black 
hole mass based on the Rbla-Lsioo relation, including nu- 
merous options for measuring the velocity width of the H/3 
emission line. These results show that the black holes masses 
derived from SED fitting have a similar distribution to that 
derived from profiles corrected for radiation pressure effects, 
except for an offset of 0.3 dex lower in both the NLS1 and 
BLS1 subsamples. The black hole mass difference between 
NLSls and BLSls from these two methods (i.e. SED fitting 
and radiation pressure corrected profiles) are both smaller 
than inferred from other mass measurements. This implies 
that compared with black hole mass estimates based only 
on the H/3 FWHM, NLSls may lie closer to the established 
M-ct* relation at the low mass end, when their black hole 
masses are corrected for radiation pressure, and when we 
use masses derived from our SED fitting. 

Finally, we form three broadband SED templates by co- 
adding SEDs in three subsamples (consist of 12 objects in 
each) to examine how the broadband SED depends on H/3 
FWHM velocity width, and by extension the Eddington ra- 
tio. The results show that there is a change in the SED shape 
as the FWHM increases, with NLSls having the largest big 
blue bump in the extreme UV region. Other important pa- 
rameters such as T2-iokeV, K2-iokeV and a ox , also change as 
the H/3 FWHM increases. The implications of correlations 
among these parameters will be discussed in our next paper. 
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Figure 12. The average SED of our sample. The panel on the left shows the averaged SED for the 12 NLSls (including two marginal 
NLSls, 2XMM 112328.0+052823 and IE 1346+26.7). The average R/3 FWHM is 1400 ± 500 km s _1 . The red area indicates a one 
standard deviation region on either side of the average spectrum. The central panel is for 12 objects with moderate line width. The 
average FWHM is 3700 ± 600 km s — 1 . The green region indicates one standard deviation. The panel on the right is the mean SED for 
the 12 broadest line objects in our sample, including the one double-peak source. The average FWHM is 9800 ± 2900 km s — 1 . We also 
show the average value of the 2-10 keV powerlaw photon index, the 2-10 keV bolometric correction, and the ct ox value with a one sigma 
error. on the Y-axis title is the luminosity distance. The unit of Y-axis is 'keV (ergs s — 1 keV - 1 )' in logarithm. The same arbitrary 
constant of 1.31 xl0 — 46 is used for rescaling each plot. 
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Figure Al. The spectral fitting results. Object order follows all other tables in this paper as increasing RA and DEC. 1. Broadband 
SED fitting plot (panel-a): X-ray data has been rebinned for each object. Green solid line is the pure accretion disc component peaking 
at optical/UV region, orange line is Comptonisation component producing soft X-ray excess below 2 keV, blue line is the hard X-ray 
Comptonisation component dominating 2-10 keV spectrum, and red is the total broadband SED model. 2. SDSS spectrum fitting plot 
(panel-b): only the fitted spectrum below 7000A is plotted. Green solid line is the best- fit underlying continuum from accretion disc. 
Orange line shows all best-fit emission lines, including the results from detailed Balmer line fitting in panel-c. Fell emission is plotted as 
light blue, while Balmer continuum being dark blue. The total best-fit model with reddening is drawn in red solid line. 3. Balmer emission 
line fitting plot (panel-c): spectral ranges containing Ha and H/3 profiles are plotted separately, with blue lines showing individual line 
components and red line showing the whole best-fit model. These are also the corresponding zoom-in plots of nearby regions of Ho and 
H/3 in panel-b. The given black hole mass is the broadband SED best-fit value, see Section [5] for detailed descriptions. 
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APPENDIX B: XMM-NEWTON AND SDSS DR7 SOURCE POSITION AND SEPERATION OF OUR 

SAMPLE 



C. Jin, M. Ward, C. Done and J. M. Gelbord 



Table Bl. XMM-Newton and SDSS DR7 source position and separation of our sample. ID: object number, the same as 
Table [T] XMM_Ra and XMM_Dec: source's right ascension and declination in the corresponding XMM-Newton observa- 
tion; XMM_PosErr: X-ray position uncertainty from XMM-Newton; SDSS_Ra and SDSS_Dec: source's right ascension and 
declination measured by SDSS; Separation: the angular separation between source's XMM-Newton and SDSS coordinates; 
Sep./XMM_PosErr: the ratio between coordinates separation and X-ray position uncertainty, showing the significance of 
coordinate separation. 
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APPENDIX C: BLACK HOLE MASSES FROM DIFFERENT METHODS 



C. Jin, M. Ward, C. Done and J. M. Gelbord 



Table CI. Black hole masses from different methods. Mjg jq: black hole mass calculated from the FWHM of H/3 in- 
termediate component in logarithm and solar mass; Mbh,bc- black hole mass calculated from the FWHM of H/3 broad 
component; VIbh IC+BC- black hole masses calculated from the FWHM of superposing H/3 intermediate component (IC) 
and broad component (BC) (i.e. narrow component subtracted), using Equation [5] Mjg CT : black hole mass calculated from 
the second momentum of the whole H/3 line profile, see subsection 16.51 for details; M-bh Fit'- the best-fit black hole masses 
in logarithm, which is constrained by M-bh,ic an d Mjj jq, but values within 0.5 lower than log(Msfr ic) were also 
ad opted in the fitting, s ee subsection 16.51 log(M. bh,rp)'- the radiation pressure corrected black hole mass using Equation 9 
in lMarconi et all l|200Sl ) with /=3.1 and logg=7.6; (*): note that Mbh,ic+bc is always within the range of Mbh,ic an d 
Mfljj ( Bc, except for UM269 whose H/3 shows double-peak profile. 
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Table Dl. Emission line parameters for the whole sample. Narrow component (NC), intermediate component (IC), broad component 
(BC) and intermediate plus broad component (I+B) are shown separately for Ha and H/3. IC and BC are both Gaussian, while NC 
may have the same profile as the whole [OIII] A5007 or only the narrowest Gaussian component in [OIII] A5007. [OIII] A5007 has two 
(or sometimes three) Gaussian components. In the case of [OIII] A5007, 'I+B' raw shows the parameters for the whole emission line 
rather than having narrow component subtracted Only one Gaussian profile is used for Hell A4686, [FeVII] A6087 and [FeX] A6374. 
Sometimes the S/N of our spectra is not high enough to resolve all these lines, or they are too week to be resolved, thus they do 
not have their line parameters measured, 'vel' means velocity of line center relative to the rest frame vacuum wavelength in ferns" . 
The velocity of 'NC is small and may come from the redshift uncertainty in Sloan's final rcdshift measurement, and thus should not 
be taken seriously. FWHMs of 'NC, 'IC and 'BC are directly from the Gaussian profile parameters. FWHM for 'I+B' is measured 
directly from the superposed model profile. The numbers are all in kms^ 1 . 'lum' and 'ew' means luminosity in Logio{ergss~ 1 ) and 
equivalent width in A. 
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401 


41.23 


12 










aic 








AO 






4117 


f, A 




— 


— 




1 1 Q 


A~\ *3 


41 9 


(l)/ 


IC 


120 


2040 


42.36 


140 


120 


2040 


41.86 


32 


52 


396 


41.59 


17 










BC 


1100 


9390 


42.05 


68 


1100 


9390 


41.80 


28 


-130 


1030 


41.51 


15 










I+B 




2190 


42.53 


210 




2310 


42.14 


59 




468 


41.85 


32 








A* 


aic 






AC\ SQ 


o.o 






AC\ ^ 


n si 


1A 


239 


40.82 


9 ^ 


AC\ 8 








IC 


-400 


5260 


42.23 


77 


-400 


5260 


41.31 


7.4 


140 


192 


41.23 


6.2 










BC 


-450 


13700 


42.29 


86 


-450 


13700 


42.20 


57 


-110 


1370 


40.87 


2.7 










I+B 




6500 


42.56 


160 




10800 


42.25 


64 




248 


41.49 


11 








5 


NC 




326 


41.76 


13 




325 


41.01 


1.6 










41.9 


41.9 




(2) / 


IC 


190 


2360 


43.15 


320 


190 


2360 


42.56 


58 


-23 


325 


41.71 


8.4 










BC 


900 


9030 


42.96 


200 


900 


9030 


42.52 


53 


-150 


794 


42.01 


17 










I+B 




2610 


43.37 


520 




2720 


42.84 


110 




474 


42.19 


25 








6 


NC 




461 


42.66 


51 




457 


41.88 


4.7 










42.4 






(1)/ 


IC 


400 


3980 


43.40 


280 


400 


3980 


42.86 


45 


-76 


387 


42.25 


12 










BC 


1400 


14400 


43.53 


370 


1400 


14400 


43.16 


90 


-120 


1020 


42.19 


10. 










I+B 




4930 


43.77 


650 




5430 


43.34 


140 




462 


42.52 


22 








7 


NC 




580 


41.88 


38 




575 


41.38 


8.0 










41.5 






(2)/ 


IC 


160 


1570 


42.10 


62 


160 


1570 


41.60 


13 


-160 


580 


40.96 


3.2 










BC 


-120 


4300 


41.99 


48 


-120 


4300 


41.70 


17 


-610 


1200 


41.27 


6.5 










I+B 




1820 


42.35 


110 




1980 


41.95 


30 




1030 


41.45 


9.7 








8* 


NC 




442 


41.75 


78 




445 


41.08 


12 


-27 


300 


41.10 


13 


41.6 


41.0 


41.0 


(1)/ 


IC 


73 


2360 


42.28 


260 


73 


2360 


41.68 


48 


-160 


663 


41.26 


19 










BC 


920 


8080 


41.97 


130 


920 


8080 


41.72 


53 


-170 


1390 


41.03 


11 










I+B 




2580 


42.45 


390 




2840 


42.00 


100 




444 


41.62 


43 








9* 


NC 




317 


41.46 


150 




319 


40.81 


34 


-10. 


297 


41.66 


240 


40.6 


40.4 




(1)/ 


IC 


200 


2360 


42.02 


530 


200 


2360 


40.99 


50 


150 


1450 


40.44 


14 










BC 


840 


7230 


41.48 


150 


840 


7230 


41.14 


72 


-41 


509 


41.10 


64 










I+B 




2500 


42.13 


680 




3030 


41.37 


120 




316 


41.79 


320 
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Table Dl. continued 



ID 






Ha 








H/9 






[OIII] 5007 




Hell 


FcVII 


FcX 






vel 


fwhm 


lum 


ew 


vel 


fwhm 


lum 


ew 


vel 


fwhm 


lum 


ew 


lum 


lum 


lum 


10 


NC 


— 


405 


42.06 


18 


— 


401 


41.13 


1.3 


— 


— 


— 


— 


41.9 


42.0 


41.2 


(iy 


IC 


180 


4060 


43.45 


430 


180 


4060 


42.87 


70 


28 


330 


41.96 


8.9 










BC 


1200 


12300 


43.12 


200 


1200 


12300 


42.83 


64 


-210 


934 


41.97 


9.2 










I+B 




4440 


43.62 


630 




4810 


43.15 


130 




407 


42.27 


18 








11 


NC 


— 


279 


42.03 


47 


— 


281 


41.38 


5.8 


— 


— 


— 


— 


41.7 


— 


— 


i-i \ f 


IC 


230 


3440 


42.89 


340 


230 


3440 


42.11 


31 


-43 


249 


42.06 


29 










BC 


730 


14100 


43.10 


550 


730 


14100 


42.66 


110 


-82 


634 


41.77 


15 










I+B 




4360 


43.31 


900 




5640 


42.77 


140 




279 


42.24 


44 








12 


NC 


— 


396 


41.55 


36 


— 


394 


40.60 


2.8 


— 


— 


— 


— 


41.6 


— 


— 


(iy 


10 


-76 


3850 


42.53 


350 


-76 


3850 


41.90 


57 


36 


1590 


40.92 


6.2 










BC 


-140 


12600 


42.11 


130 


-140 


12600 


41.77 


43 


-23 


372 


41.33 


16 










I+B 




4130 


42.67 


480 




4390 


42.14 


100 




395 


41.47 


22 








13 


NC 


— 


396 


42.50 


10. 


— 


395 


41.94 


1.6 


— 


— 


— 


— 


42.8 


42.3 


— 


n \ f 

(iy 


10 


580 


6160 


43.53 


110 


580 


6160 


42.85 


13 


-44 


304 


42.16 


2.8 










BC 


700 


17200 


43.62 


140 


700 


17200 


43.46 


52 


-250 


1190 


42.37 


4.5 










I+B 




7730 


43.88 


240 




10800 


43.55 


65 




395 


42.58 


7.3 








14 


NC 


— 


298 


41.97 


23 


— 


300 


41.48 


5.8 


— 


— 


— 


— 


42.0 


41.9 


— 


n \ f 

(iy 


10 


300 


5580 


43.31 


520 


300 


5580 


42.59 


74 


6.0 


893 


42.10 


25 










BC 


1900 


16300 


42.81 


160 


1900 


16300 


42.70 


96 


-49 


266 


42.31 


41 










I+B 




5940 


43.43 


680 




7060 


42.95 


170 




297 


42.52 


65 








15 


NC 


— 


342 


41.03 


30 


— 


345 


40.43 


7.4 


— 


— 


— 


— 


40.7 


40.0 


40.4 


(2y 


10 


58 


918 


41.22 


46 


58 


918 


40.73 


15 


53 


299 


40.84 


19 










BC 


-20 


4400 


40.97 


26 


-20 


4400 


40.48 


8.3 


-320 


1000 


40.74 


15 










I+B 




994 


41.42 


73 




987 


40.92 


23 




340 


41.09 


34 








16 


NC 


— 


279 


41.62 


9.9 


— 


281 


40.98 


1.5 


— 


— 


— 




42.1 


41.4 


— 




10 


37 


2810 


43. 14 


330 


37 


2810 


42.47 


46 


-55 


968 


41.53 


0.4 










BC 


1300 


8860 


42.85 


170 


1300 


8860 


42.61 


63 


-76 


280 


42.00 


16 










I+B 




3080 


43.32 


500 




3560 


42.84 


110 




297 


42.13 


22 








17* 


NC 


— 


608 


42.00 


37 


— 


607 


41.08 


2.9 


-28 


617 


41.84 


17 


41.7 


41.8 


41.7 


/o\ f 

(3) J 


IC 


160 


1930 


42.91 


300 


160 


1930 


42.35 


53 


140 


235 


41.14 


3.4 










BC 


430 


7450 


42.67 


170 


430 


7450 


42.34 


51 


-530 


1600 


41.57 


9.2 










I+B 




2120 


43.10 


470 




2250 


42.64 


100 




607 


42.08 


30 








18 


NC 


— 


400 


42.23 


32 


— 


401 


41.53 


3.6 


— 


— 


— 


— 


42.1 


— 


41.9 


(1)J 


10 


220 


1810 


43.02 


200 


220 


1810 


42.44 


29 


-60 


260 


41.28 


2.1 










BC 


360 


7220 


42.93 


160 


360 


7220 


42.66 


48 


-210 


701 


41.61 


4.6 










I+B 




2060 


43.28 


360 




2310 


42.86 


76 




401 


41.78 


6.7 








19 


NC 




188 


40.90 


13 




191 


40.25 


2.3 




— 


— 


— 


41.0 


40.4 


40.5 


(2) J 


10 


A r> 
ID 


loOU 


1 "1 77 

41. 1 l 


96 


A a 

4b 


1500 


41.16 


i n 
19 


to 


186 


40.93 


11 










BC 


-54 


5010 


41.66 


74 


-54 


5010 


41.40 


33 


-100 


445 


40.80 


8.4 










I+B 




1730 


42.02 


170 




2000 


41.60 


51 




223 


41.17 


20 








20 


NC 


— 


232 


40.64 


14 


— 


230 


40.03 


3.1 


— 






— 


41.2 


— 


39.8 


(2)/ 


IC 


29 


578 


41.12 


40 


29 


578 


40.51 


9.5 


43 


234 


40.29 


5.8 










BC 


120 


1970 


41.12 


41 


120 


1970 


40.76 


17 


-230 


460 


39.80 


1.9 










I+B 




686 


41.42 


82 




774 


40.95 


27 




254 


40.41 


7.6 








21 


NC 




354 


42.05 


49 




357 


41.53 


11 












41.5 




(2)/ 


IC 


-250 


3990 


43.03 


470 


-250 


3990 


42.14 


44 


110 


357 


42.29 


64 










BC 


730 


13900 


42.93 


380 


730 


13900 


42.58 


120 


-200 


1160 


42.27 


63 










I+B 




4570 


43.29 


840 




6090 


42.71 


170 




419 


42.58 


130 








22 


NC 




317 


40.86 


6.5 




319 


40.08 


0.86 










41.5 






(2)/ 


IC 


750 


4700 


42.28 


170 


750 


4700 


41.59 


28 


-19 


316 


40.92 


6.0 










BC 


59 


12600 


42.18 


130 


59 


12600 


41.98 


68 


-150 


764 


40.90 


5.8 










I+B 




5510 


42.53 


300 




7050 


42.13 


96 




389 


41.21 


12 









A Spectral Study of Unobscured Type 1 AGN - I 45 



Table Dl. continued 



ID 






Ha 








H/9 








[OIII] 5007 




Hell 


FcVII 


FcX 






vel 


fwhm 


lum 


ew 


vel 


fwhm 


lum 


ew 


vel 


fwhm 


lum 


ew 


lum 


lum 


lum 


23 


NC 


— 


340 


41.56 


39 


— 


344 


40.79 


5.0 


— 


— 


— 


— 


41.0 


— 


41.0 


(iy 


IC 


130 


1700 


42.15 


150 


130 


1700 


41.66 


37 


31 


256 


40.59 


3.2 










BC 


340 


5970 


41.92 


89 


340 


5970 


41.60 


33 


-79 


930 


40.81 


5.4 










I+B 




1890 


42.35 


240 




1980 


41.93 


70 




340 


41.01 


8.7 








24 


NC 


— 


576 


41.37 


3.9 


— 


575 


41.04 


1.6 


— 


— 


— 


— 


42.3 


42.2 


— 


l-l \ f 

(iy 


IC 


150 


4700 


42.74 


91 


150 


4700 


41.59 


5.5 


44 


327 


41.00 


1.4 










BC 


1400 


15400 


43.19 


260 


1400 


15400 


42.91 


120 


-110 


1090 


41.52 


4.8 










I+B 




7280 


43.32 


350 




13900 


42.93 


120 




577 


41.64 


6.2 








25 


NC 


— 


400 


42.49 


47 


— 


395 


41.79 


5.6 


— 


— 


— 


— 


42.1 


— 


— 


\ f 

(iy 


10 


-74 


3820 


43.36 


350 


-74 


3820 


42.72 


47 


29 


365 


42.51 


31 










BC 


-87 


12700 


43.22 


250 


-87 


12700 


42.92 


75 


-99 


943 


42.13 


13 










I+B 




4350 


43.59 


590 




4980 


43.13 


120 




395 


42.66 


43 








26 


NC 


— 


543 


42.16 


88 


— 


544 


41.63 


16 


— 


— 


— 


— 


41.5 


— 


— 


( 2 ) J 


10 


44 


1540 


42.53 


210 


44 


1540 


42.05 


43 


-40 


545 


41.71 


21 










BC 


490 


6770 


42.42 


160 


490 


6770 


41.92 


32 


-330 


1630 


41.21 


6.6 










I+B 




1730 


42.78 


370 




1720 


42.29 


76 




577 


41.83 


27 








27* 


NC 


— 


255 


40.98 


7.6 


— 


256 


40.48 


1.8 


49 


254 


41.53 


20 


41.5 


40.9 


40.9 


(2)1 


10 


-150 


2500 


42.40 


200 


-150 


2500 


41.64 


26 


-130 


504 


41.15 


8.6 










BC 


780 


7780 


42.43 


220 


780 


7780 


42.17 


86 


-350 


1240 


41.00 


6.1 










I+B 




3030 


42.72 


420 




4310 


42.28 


110 




291 


41.76 


35 








28 


NC 


— 


368 


42.31 


55 


— 


369 


41.74 


11 


— 


— 


— 


— 


42.7 


— 


— 


(iy 


10 





3260 


42.99 


260 





3260 


42.32 


43 


-9.0 


356 


42.62 


88 










BC 


-690 


10900 


42.75 


150 


-690 


10900 


42.52 


68 


-21 


1530 


41.94 


18 










I+B 




3610 


43.19 


420 




4240 


42.74 


110 




365 


42.71 


110 








29 


NC 


— 


497 


42.03 


34 


— 


494 


41.40 


5.6 


— 


— 


— 


— 


41.8 


41.5 


41.0 


(1)' 


10 


12 


2810 


42.87 


240 


12 


2810 


42.30 


44 


39 


385 


41.69 


11 










BC 


810 


9770 


42.71 


170 


810 


9770 


42.46 


64 


-160 


789 


41.75 


13 










I+B 




3170 


43.10 


400 




3560 


42.69 


110 




498 


42.02 


24 








30 


NC 


— 


386 


41.31 


39 


— 


388 


40.66 


5.1 


— 


— 


— 


— 


41.1 


39.8 


40.3 


ci ~\ f 

(iy 


IC 


110 


808 


41.43 


52 


110 


808 


41.09 


14 


-19 


336 


40.97 


11 










BC 


110 


3040 


41.41 


49 


110 


3040 


41.10 


14 


-210 


703 


40.74 


6.5 










I+B 




943 


41.72 


100 




953 


41.39 


28 




389 


41.17 


17 








31 


NC 


— 


382 


41.61 


8.0 


— 


382 


40.54 


0.51 


— 


— 


— 


— 


41.2 


— 


— 


(3) J 


10 


610 


5730 


43.14 


270 


610 


5730 


42.59 


57 


20 


333 


41.29 


3.0 










BC 


850 


21500 


42.75 


110 


850 


21500 


42.62 


61 


-460 


814 


41.17 


2.3 










I+B 




6140 


43.29 


380 




6810 


42.90 


120 




377 


41.54 


5.2 








32 


NC 


— 


425 


41.24 


23 


— 


424 


40.61 


5.3 


— 


— 


— 




41.1 


41.0 


40.3 


l^ \ f 
(1)1 


10 


140 


2550 


4Z. 1 1 


200 


140 


2550 


41.48 


39 


-58 


388 


41.51 












BC 


1000 


8070 


41.68 


63 


1000 


8070 


41.52 


43 


-130 


1170 


41.09 


16 










I+B 




2710 


42.30 


260 




3100 


41.80 


82 




419 


41.65 


58 








33 


NC 


— 


423 


41.68 


16 


— 


426 


41.03 


2.7 


— 






— 


41.6 


— 


— 


(1)/ 


IC 


17 


4700 


42.92 


290 


17 


4700 


42.37 


59 


100 


412 


41.99 


25 










BC 


23 


16700 


42.66 


160 


23 


16700 


42.44 


68 


-610 


1370 


41.42 


6.8 










I+B 




5170 


43.11 


450 




5690 


42.71 


130 




425 


42.09 


32 








34 


NC 




838 


42.57 


44 




836 


41.70 


4.0 










41.8 






(3)' 


IC 


260 


2470 


43.06 


130 


260 


2470 


42.52 


26 


-150 


1440 


42.48 


25 










BC 


1200 


6570 


42.97 


110 


1200 


6570 


42.72 


41 


-240 


634 


42.30 


16 










I+B 




2910 


43.32 


250 




3310 


42.93 


67 




829 


42.70 


41 








35 


NC 




572 


42.67 


130 




576 


42.01 


22 










42.0 


41.8 


42.1 


(1)/ 


IC 


270 


2380 


43.31 


560 


270 


2380 


42.65 


96 


-120 


1250 


42.28 


43 










BC 


1300 


5880 


43.01 


280 


1300 


5880 


42.54 


76 


-140 


512 


42.55 


79 










I+B 




2640 


43.49 


840 




2790 


42.90 


170 




571 


42.74 


120 
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Table Dl. continued 



ID 






Ha 








H/9 








[OIII] 5007 




Hell 


FeVII 


FcX 






vel 


fwhm 


lum 


ew 


vel 


fwhm 


lum 


ew 


vel 


fwhm 


lum 


ew 


lum 


lum 


lum 


36* 


NC 


— 


274 


41.15 


36 


— 


271 


40.32 


5.1 


54 


175 


40.54 


8.6 


41.1 


40.2 


39.6 


(1)1 


IC 


210 


1540 


41.29 


49 


210 


1540 


40.69 


12 


-62 


448 


40.62 


10. 










BC 


-160 


4450 


40.98 


24 


-160 


4450 


40.75 


14 


-290 


1150 


40.42 


6.5 










I+B 




1690 


41.47 


73 




1890 


41.02 


26 




248 


41.01 


25 








37 


NC 


— 


567 


41.92 


13 


— 


566 


41.17 


2.2 


— 


— 


— 


— 


42.5 


41.7 


41.0 


n \ f 

(iy 


IC 


220 


3440 


43.17 


230 


220 


3440 


42.68 


72 


72 


294 


41.37 


3.7 










BC 


870 


10000 


42.88 


120 


870 


10000 


42.55 


53 


-170 


824 


41.79 


9.8 










I+B 




3790 


43.35 


350 




3960 


42.92 


130 




558 


41.93 


13 








38* 


NC 


— 


326 


41.61 


140 


— 


325 


41.05 


38 


38 


635 


41.22 


56 


41.1 


— 


— 


n \ f 

(iy 


iC 


-100 


4830 


42.23 


580 


-100 


4830 


41.25 


60 


60 


286 


41.51 


110 










BC 


-390 


12200 


41.83 


230 


-390 


12200 


41.51 


110 


-140 


1280 


40.77 


20 










I+B 




5250 


42.38 


800 




6630 


41.70 


170 




328 


41.74 


190 








39 


NC 


— 


212 


41.00 


32 


— 


217 


40.44 


6.9 


— 


— 


— 


— 


40.8 


40.0 


40.1 


l-t \ f 

(iy 


10 


29 


829 


41.47 


93 


29 


829 


40.87 


19 


33 


193 


41.06 


29 










BC 


50 


3000 


41.25 


57 


50 


3000 


40.90 


20 


-20 


427 


40.73 


14 










I+B 




925 


41.67 


150 




990 


41.19 


38 




216 


41.22 


43 








40* 


NC 


— 


446 


41.34 


8.0 


— 


451 


40.45 


0.88 


-130 


465 


40.73 


1.7 


— 


— 


41.0 


[3y 


IC 


210 


1930 


42.64 


160 


210 


1930 


41.84 


22 


130 


123 


40.59 


1.2 










BC 


76 


6350 


42.47 


110 


76 


6350 


42.21 


51 


-160 


1590 


41.27 


6.0 










I+B 




2180 


42.86 


270 




2790 


42.37 


73 




450 


41.45 


8.9 








41 


NC 


— 


400 


42.10 


23 


— 


401 


41.47 


3.5 


— 


— 


— 


— 


42.3 


— 


— 


(1\ f 

(2) / 


iO 


170 


2300 


43.38 


440 


170 


2300 


42.85 


84 


-52 


401 


42.32 


25 










BC 


890 


8300 


43.09 


220 


890 


8300 


42.73 


62 


-260 


908 


42.25 


22 










I+B 




2510 


43.56 


660 




2610 


43.10 


150 




486 


42.59 


47 








42 


NC 


— 


307 


40.74 


14 


— 


306 


39.90 


1.6 


— 


— 


— 


— 


41.1 


40.6 


39.9 




iO 


120 


3760 


41.89 


190 


120 


3760 


41.27 


38 


26 


721 


40.43 


5.6 










BC 


-300 


9930 


41.66 


110 


-300 


9930 


41.44 


56 


-3.0 


287 


41.02 


22 










I+B 




4240 


42.09 


300 




4920 


41.67 


94 




303 


41.12 


27 








43 


NC 


— 


414 


42.10 


29 


— 


419 


41.60 


6.5 


— 


— 


— 


— 


42.0 


— 


41.0 


(iy 


IC 


12 


3550 


43.23 


380 


12 


3550 


42.51 


53 


79 


324 


42.07 


20 










BC 


1400 


10900 


42.95 


200 


1400 


10900 


42.67 


75 


-310 


993 


42.22 


28 










I+B 




3920 


43.41 


590 




4550 


42.90 


130 




419 


42.45 


47 








44 


NC 




265 


41.86 


26 




262 


41.34 


5.1 


— 


— 


— 


— 


41.8 


40.3 


40.8 


(iy 


iO 


(6 


852 


42.41 


91 


1 3 


852 


41.80 


15 


110 


220 


41.67 


11 










BC 


97 


4190 


42.32 


74 


97 


4190 


42.04 


26 


-98 


852 


41.72 


13 










I+B 




952 


42.67 


170 




1070 


42.24 


40 




260 


41.99 


24 








45 


NC 


— 


345 


42.82 


33 


— 


344 


42.22 


5.2 
















(iy 


iO 


-870 


6210 


43.81 


330 


-870 


6210 


43.01 


32 


15 


267 


42.79 


20 










BC 


1400 


17200 


43.89 


390 


1400 


17200 


43.56 


110 


-62 


880 


42.93 


28 










I+B 




7740 


44.16 


720 




10900 


43.66 


140 




346 


43.17 


48 








46* 


NC 


— 


514 


42.22 


51 


— 


516 


41.66 


13 


-11 


432 


42.32 


63 


41.4 


41.9 


— 


(3) / 


IC 


3000 


6960 


42.54 


110 


3000 


6960 


41.99 


29 


150 


223 


41.66 


14 










BC 


-1400 


8410 


42.93 


260 


-1400 


8410 


42.09 


36 


-150 


1010 


42.30 


60 










I+B 




9970 


43.08 


370 




9930 


42.34 


64 




534 


42.66 


140 








47 


NC 




391 


41.58 


27 




391 


40.93 


5.7 










41.7 






(1)/ 


IC 


160 


3550 


42.70 


360 


160 


3550 


41.94 


58 


50 


1080 


41.19 


11 










BC 


670 


7400 


42.02 


74 


670 


7400 


41.78 


40 


-24 


355 


41.70 


34 










I+B 




3730 


42.78 


430 




4100 


42.17 


98 




383 


41.82 


45 








48 


NC 




302 


41.27 


51 




306 


40.67 


8.7 










40.4 




40.2 


(2)/ 


IC 


110 


1040 


41.53 


92 


110 


1040 


41.03 


20 


28 


306 


40.37 


4.6 










BC 


33 


4300 


41.38 


66 


33 


4300 


41.00 


19 


-130 


719 


40.30 


3.9 










I+B 




1160 


41.77 


160 




1190 


41.32 


39 




365 


40.64 


8.4 
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Table Dl. continued 



ID 






Ha 








H/3 








[OIII] 5007 




Hell 


FeVII 


FeX 






vel 


fwhm 


lum 


ew 


vel 


fwhm 


lum 


ew 


vel 


fwhm 


lum 


ew 


lum 


lum 


lum 




jv u 




Zo4 


41. Z i 


1 

iy 




OO 1 
iOl 


4U.D0 


i a 
o . o 










41.0 




4U.0 


(1)/ 


IC 


16 


1010 


41.65 


15 


16 


1010 


41.05 


9.4 


11 


268 


40.98 


8.1 










BC 


160 


4040 


41.49 


31 


160 


4040 


41.33 


18 


-260 


578 


40.37 


2.0 










I+B 




1120 


41.88 


76 




1340 


41.51 


27 




285 


41.07 


10. 








50* 


NC 




274 


42.84 


36 




275 


42.23 


4.5 


-73 


402 


41.87 


2.1 


42.3 


42.2 


41.9 


(1)/ 


IC 


120 


1880 


43.60 


210 


120 


1880 


43.06 


30 


130 


224 


42.34 


6.2 










BC 


91 


8960 


43.50 


170 


91 


8960 


43.12 


35 


-470 


1770 


42.42 


7.4 










I+B 




2090 


43.86 


370 




2200 


43.39 


65 




266 


42.74 


16 








51* 


NC 




451 


41.73 


65 




445 


41.20 


16 


11 


364 


42.03 


110 




41.2 




(1)/ 


IC 


-1300 


6490 


42.44 


330 


-1300 


6490 


41.59 


39 


130 


1150 


41.57 


39 










BC 


1400 


17200 


42.53 


410 


1400 


17200 


42.10 


130 


340 


247 


41.23 


18 










I+B 




8170 


42.79 


750 




11100 


42.21 


170 




450 


42.21 


170 









f : The final fitting method used for each object, see subsection 13.11 for detailed description of each fitting methods. 
* : Three gaussian profiles are used for these objects, in this case (I+B) means the total of all three components. 

d : UM 269, the only object in our sample showing double-peak feature in Balmer lines. Two gaussian profiles are used for fitting the 
two peaks, thus the velocity shift of each component related to the line centre is huge. 



